A FREE BOUNDARY PROBLEM INSPIRED BY A CONJECTURE OF DE 

GIORGI 



Abstract. We study global monotone solutions of the free boundary problem that arises from 
minimizing the energy functional I(u) = J |Vu| 2 + V(u), where V(u) is the characteristic function 
of the interval (—1, 1). This functional is a close relative of the scalar Ginzburg-Landau functional 
J(«) = J | V«| 2 + W(u), where W(u) = (1 — u 2 ) 2 /2 is a standard double-well potential. According 
to a famous conjecture of De Giorgi, global critical points of J that are bounded and monotone 
in one direction have level sets that are hyperplanes, at least up to dimension 8. Recently, Del 
Pino, Kowalczyk and Wei gave an intricate fixed-point-argument construction of a counterexample 
in dimension 9, whose level sets "follow" the entire minimal non-planar graph, built by Bombieri, 
De Giorgi and Giusti (BdGG). In this paper we turn to the free boundary variant of the problem 
and we construct the analogous example; the advantage here is that of geometric transparency 
as the interphase {\u\ < 1} will be contained within a unit-width band around the BdGG graph. 
Furthermore, we avoid the technicalities of Del Pino, Kowalczyk and Wei's fixed-point argument 
by using barriers only. 
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1. Introduction. 
In this paper we study the following free boundary problem: 

Au — in n iD :— {x e f2 : |u(a;)| < 1} 
u = ±l in 0\n in (1.1) 
|Vu|=l on dn in n Q = F + (u) U F~(u) 

where f2 C K™ is a domain and the free boundary of u consists of two pieces 

F+{u) := d{u = 1} n n F~(u) := <9{m = -1} n ft. 

In particular, we are interested in global solutions (ft = R n ), which are monotonically increasing in 
the last coordinate x n . We pose the following question: 

Problem *. Let n = 9. Does i/iere exist a global solution to (1.1), monotonically increasing in xg, 
such that its level sets are not hyperplanes? 

The question above should be read in the context of the prominent De Giorgi's conjecture con- 
cerning global solutions of the Allen-Cahn equation 

Au = (l-u> inM n , (1.2) 

namely: 

Conjecture (De Giorgi [ ]). If u 6 C 2 (R n ) is a global solution of (1.2) such that \u\ < 1 and 
d Xn u > 0, then the level sets {u = A} are hyperplanes, at least for dimensions n < 8. 

The common nature of the PDE's (1.1) and (1.2) is rooted in the fact that they arise as Euler- 
Lagrange equations for the closely related energy functionals I and J, respectively 

I(u,Q)= I |Vu| 2 + V(u) for u:fi->.[-l,l] (1.3) 
Jn 

J(u,Q) = [ |Vu| 2 + W(«), (1.4) 
Jn 

where V(u) := l(_i t i\(u) is a singular version of the standard double-well potential W{u) := - — ■ 
De Giorgi's conjecture has been motivated by a fascinating and deep connection between the theory 
of semilinear elliptic PDE and the theory of minimal surfaces. The connection was first rigorously 
stated through the notion of F-convergence in the work of Modica [16]. Assuming that u minimizes 
J in a large ball -Bi/ e , u e (x) = u{x/e) minimizes the rescaled energy 

J e (v,Bi) = e / \\7v\ 2 + - I W{v) 

J Bi e JBt 

in the unit ball. What Modica proved was that as e — ¥ 0, a subsequence of minimizers it £fc of J £fc (•, B\) 
of uniformly bounded energy converges to 

Ue k -+ 1e ~ IbiVe in iL^i), 

where E has a perimeter minimizing boundary in Bi, i.e. dE is a minimal hypersurface. The 
convergence is in fact stronger: Caffarelli and Cordoba ([ ], [ ]) later showed that the level sets 
{u e — A} for — 1 < A < 1 converge uniformly on compacts to the minimal hypersurface dE. Intuitively 
speaking therefore, the level sets of a global minimizer of J look like a minimal hypersurface at large 
scales. The analogous statements can, of course, be made for global minimizers of /. 
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The monotonicity assumption d Xn u > in De Giorgi's conjecture implies that u is not only a 
stable critical point for J, but that u is, in fact, a global minimizer of J in a certain sense (see [1]). 
Under the natural assumption 

lim u(x', x n ) = ±1. (1-5) 

x n — >-±oo 

the level sets of u are also graphs over K ra_1 in the e„-direction. Therefore, the preceding discussion 
combined with Bernstein's theorem, which states that an entire minimal graph in IR™ -1 x R is a 
hyperplane for n < 8 (cf. Simons [ ]), is what gives plausibility to De Giorgi's conjecture. Moreover, 
the existence of a non-planar minimal graph in dimension n — 9, constructed by Bombieri, De Giorgi 
and Giusti in [5], strongly suggests that the conjecture is likely to be false for n > 9. 

There has been a lot of recent work which has almost completely resolved De Giorgi's conjecture. 
It was fully established in dimensions n = 2 by Ghoussoub and Gui [13] and n = 3 by Ambrosio 
and Cabre [4], while Savin [18] managed to prove it for dimensions 2 < n < 8 under the additional 
assumption (1.5). Savin's approach has a broader scope and applies to monotone global minimizers 
of the functional /, as well (see [19]). 

Recently Del Pino, Kowalczyk and Wei [1 2] successfully constructed a counterexample in dimension 
n = 9. Roughly speaking, their strategy is based upon the derivation of a sufficiently good ansatz 
whose level sets "follow" the Bombieri-De Giorgi-Giusti (BdGG) minimal graph. This allows them 
to carry out an intricate fixed point argument. 

It was a desire to gain a better understanding of precisely what geometric ingredients are respon- 
sible for the existence of this important counterexample that led us to formulate and resolve in the 
affirmative the alternative Problem *. 

Theorem 1.1. There exists a solution u : K 9 — > M. of (1.1) which is monotonically increasing in xg 
and whose free boundary F(u) = F + (u) U F~(u) consists of two non-planar smooth graphs. 

The study of the free boundary problem has an obvious geometric advantage. In this setting, the 
interphase {\u\ < 1} will be contained within a unit-width band around the BdGG minimal graph, 
so that ones does not have to worry about capturing a non-trivial behaviour of the solution away 
from the band. We will use the method of barriers which is elementary in nature and allows for a 
transparent and relatively precise description of the solution (we will be able to trap the solution quite 
tightly between the two barriers). This way we avoid using fixed point arguments which are arguably 
the main culprit for the level of technical complexity of the construction by Del Pino, Kowalczyk and 
Wei. 

To construct a solution to (1.1) once we are in possession of a supersolution lying above a subso- 
lution (we will define these notions shortly), we adopt the strategy developed by De Silva [10] in her 
study of global free boundary graphs that arise from monotone solutions to the classical one-phase 
free boundary problem: 

Au = in Ct p (u) := {x G Q : u{x) > 0} , , 

|Vu| = 1 on F p {u) :=Qr\dn p (u). ^ ' ' 

Namely, a global solution to (1.6) is constructed as the limit of a sequence of local minimizers of the 
one-phase energy functional: 

7 (tt,n) - f \Vu\ 2 + 1 {U>0} (1.7) 
Jn 

constrained to lie between a fixed strict subsolution and a fixed strict supersolution to (1.6). The 
strictness condition ensures that the free boundary of each minimizer doesn't touch the free bound- 
aries of the barriers. Following the classical ideas of Alt, Cafarelli and Friedman ([2], [3]), De Silva 
shows that u is a global energy minimizing viscosity solution, which is locally Lipschitz continuous 
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and has non-degenerate growth along its free boundary; moreover, if one assumes that the two bar- 
riers are monotonically increasing in x n , the global solution can also be chosen to be monotonically 
increasing in x n after a rearrangement. The harder part is the regularity theory: De Silva's key 
observation is that the positive phase of the minimizer is locally an NTA (non-tangentially acces- 
sible) domain ([14]) which allows the application of the powerful boundary Harnack principle. By 
comparing the solution with a vertical translate, she rules out the possibility that the free boundary 
contains any vertical segments, so that it is a graph, and then she shows that the graph is, in fact, 
continuous. A more sophisticated comparison argument by De Silva and Jerison [11] establishes a 
Lipschitz bound on the free boundary graph. Hence, by the classical result of Caffarelli [6], the free 
boundary is locally a C 1,a graph, so that the global minimizer is indeed a classical solution to (1.6). 

Obviously, the functionals I and Jo are close relatives: if u minimizes /(-,f2) and D C 57 is a 
(nice enough) subdomain, such that D n {u = 1} = 0, then u+ 1 minimizes Iq(-,D); similarly, if 
Dn{w = —1} = 0, 1 — u minimizes Zo(-, D). So, after we construct a global minimizer u to /, we will 
be in a position to directly apply the regularity theory for the one-phase energy minimizers from the 
discussion above to the free boundary of u. 

Let us now give a brief outline of the arguments and the structure of our paper. 



First, let us recall the definition of a classical super/subsolution to the one-phase problem (1.6) 
(see for example [6]). 

Definition 2.1. A classical super solution (resp. subsolution) to (1.6) is a non-negative continuous 
function w in Q such that 

• we c 2 {n p {w)) ; 

• Aw < (resp. Aw > 0) in Cl p (w); 

• The free boundary F p (w) is a C 2 surface and 



If the inequality above is strict, we call w a strict super (resp. sub) solution. 

The appropriate notion of a classical super/subsolution to our free boundary problem (1.1) is, 
therefore, the following: 

Definition 2.2. A classical supersolution (resp. subsolution) to (1.1) is a non-negative continuous 
function w in 51 such that 

• w£ C 2 (n m (w)); 

• Aw < (resp. Aw > 0) in fl in (w); 

• The free boundary F(w) = F + (w) U F~(w) consists of two C 2 surfaces and 



If the inequalities above are strict, we call w a strict super (resp. sub) solution to (1.1). 

As mentioned in the introduction, the driving intuition is that the level surfaces of a solution to 
(1.1) should follow the shape of the BdGG entire minimal graph T = {(x', xg) € K 8 x K : xg — F(x')}. 
The function F : M 8 —> R satisfies the minimal surface equation H[F] = where H[-j is the mean 
curvature operator (MCO) 



2. Outline of strategy. 



< |Vw| < 1 (resp. |Viu| > 1) on F p (w). 



< |Vw| < 1 (resp. \Vw\ > 1) on F~(w), 
\Vw\ > 1 (resp. \Vw\ < 1) on F + (w). 
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Note that there is a whole one-parameter family of such entire minimal non-planar graphs, obtained 
by rescaling F: 

r Q := a" x r = {x 9 = F a (x')} 

where a > and F a (x') := a~ 1 F(ax'). "Following the shape" should be interpreted in the sense 
that we would like the solution u and thus the trapping super /subsolution W and V to behave 
asymptotically at infinity like the signed distance to r Q for some a > 0: 

V{x) < u(x) < W(x) « signed dist(a;,r„) 

within their interphases (for the definition of the signed distance: we take the sign to be positive if 
the point x e R 9 lies "above the graph", i.e. if xg > F a (x'), and negative otherwise). This suggests 
that the coordinates 

R 9 3 x -> (y, z) £ r Q x K, x = y + zv a (y), (2.1) 

where v a (y) is the unit normal to r a at y with v a {y) • eg > 0, will be particularly well-suited to the 
problem. We will later show in Lemma 3.1 that the coordinates 

R 9 3 x -> (y, z) G r x x K, x = y + zv 1 (y), (2.2) 

are well-defined in a thin band around T = Tj 

B r (d) ={i€l 9 : dist(x,r) < d} 

for d > small enough. By taking blow-ups of space x — > we see the coordinates (2.1) with 

respect to the blow-up T a will be well defined in the band Br a {ct~ 1 d). Thus we can ensure that the 
coordinate system (2.1) is well defined in a unit-size band Br a = 2?r Q (2) for all a > small enough. 
The trick of scaling will prove quite useful in what comes later, as well. Its effect on the geometry in 
the unit-width band is described in §3, Lemma 3.4. 

Not surprisingly, we look for a supersolution/subsolution that is a polynomial in z: 

m 

w(y,z) = Y, h k(y> k , 

k=0 

where hf ~ 1 to main order at infinity, whereas all the other coefficient decay appropriately to 0. As 
it turns out, m = 5 suffices. 

The Euclidean Laplacian is given by the following key formula: 

A = A raiz) + d*-H ra{z) (y)d z , 

where As denotes the Laplace-Beltrami operator on a surface S, 

F a (z) = {y + zv a {y) : ye T} 

is a level set for the signed distance to T Q and Hj;r z -\(y) is its mean curvature at y + zv(y). Note 
that if kf (y) denote the principal curvatures of T a at y 

Expanding (2.3) in 2, we formally have 

00 

Hr a(z )(y) = (H Ul =Q) + Y J ^ 1 H Ua , (2.4) 

1=2 

where Hi a — Yl--,(kf) 1 is the sum of the l-ih powers of the principal curvatures of T a . It turns out 
that the principal curvatures 

kf(y) = 0(a(l + alx'l)- 1 ), y = {x',F a {x')) G F Q 
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meaning in particular that the series (2.4) converges rapidly for (y, z) 6 Br a and a small enough. 
However, a more refined understanding of the asymptotics of the quantities H^ := H^i will be 
needed. 

Del Pino, Kowalczyk and Wei faced the exact same issue in [12]. For the purpose, they introduce 
the model graph L^, which has an explicit coordinate description and which matches T very well 
at infinity. This allows one to approximate geometric data and geometric operators on T with their 
counterparts on T^: for example, the second fundamental form, the quantities Hi, the intrinsic 
gradient and the Laplace-Beltrami operator. We will briefly present their results concerning the 
geometry of T and the closeness between T and Too, and use their framework to prove additional 
relevant estimates in Section 3. 

Choosing the coefficients h%{y) so that w meets the supersolution conditions is the subject of 
Section 4. In fact, most of the choices will be imposed on us (see Remark 4.1): they will be given 
in terms of geometric quantities like Hi, a and their covariant derivatives. The upshot is that (see 
Lemmas 4.1 and 4.3) 

Aw = Jr hn + h'j* — z 2 H-\ a + lower order terms in Br , 

a 2 „, (25) 

|Vtu| = 1 ± h£ + lower order terms on {w = ±1}, 

where 

J Fa = A Fa + \A a \ 2 

is the Jacobi operator on T a and h'g = ZW; — \A a \ 2 hQ. Thus, by varying Hq we can satisfy the 
superharmonicity condition and by varying h'^ - the gradient condition on the free boundary. So, 
h'2 needs to be positive and we will also require that hg > 0. That way, we only have to flip the 
signs of the coefficients and in the ansatz (and leave the remaining ones unchanged) in order 
to produce a subsolution ansatz that automatically lies underneath the supersolution. 

So, we want the function /iq to be a positive supersolution for Jr a that satisfies an appropriate 
differential inequality. It turns out that Jr admits nonnegative supersolutions h of the following 
types: 

• Type 1 is such that Jrh can absorb terms that decay like r~ k for k > 4. See Proposition 4.1. 
This is useful when dealing with the lower order terms in (2.5). 

• Type 2 can take care of the |i?3|-term which is globally 0(r~ 3 ) but has the important 
additional property that it vanishes on the Simons cone S = {(u,v) e R 4 x I 4 : \u\ = \v\}. 

The Type 1 supersolution is readily provided by Del Pino, Kowalczyk and Wei's paper [12, Proposi- 
tion 4.2(b)]. We build the Type 2 supersolution ourselves in Section 4.2 (asymptotically in Lemma 
4.6 and globally in Proposition 4.2) and the construction involves a very delicate patching of two 
supersolutions in a region of the graph over the Simons cone. The ingredients are contained in the 
analysis of the linearized mean curvature operator H'lF^] around Too carried out in [12, §7], Lemma 
7.2 and 7.3; for the reader's convenience we state these results in Appendix A. 3. The operators 
Jr^ := Ar^ + |^4oo| 2 and H'[Foo] are closely related (see (A.l)) and in turn Jr^ is asymptotically 
close to Jr (see (3.24)). This allows one to first build a (weak) supersolution away from the origin, 
which can then be upgraded to a global smooth supersolution via elliptic theory. 

Having these two types of barriers for Jr we will be able to satisfy the free boundary supersolution 
conditions far away from the origin. In order to satisfy them globally, we employ the trick of scaling 
by a. That way we can also ensure that both the supersolution and the subsolution are monotonically 
increasing in xg for all small enough a > 0. 

Once we have obtained the monotone sub-super solution pair V < W we proceed to construct 
the solution u to (1.1) as a global minimizer u of the functional /, constrained to lie in-between 
V < u < W. This is the subject of Section 5. 
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3. The Bombieri-De Giorgi-Giusti graph r and its approximation. 

In this section we will describe several results concerning the asymptotic geometry of the entire 
minimal graph in 9 dimensions, constructed by Bombieri, De Giorgi and Giusti in [•")]. Some of them 
have been covered by the analysis of the graph T, carried out by Del Pino, Kowalczyk and Wei 
in [12] (Lemmas 3.1, 3.2, 3.4 and 3.5 below). We will state those in a form suitable for our later 
computations. Furthermore, we will establish the important estimates for the covariant derivatives 
of Hi in Lemma 3.3. 

Let us first set notation. Recall, we denote the entire minimal graph by 

T = {(V, xg) e R 8 x I : i 9 = F(x')} C I 8 x I, 
where F : R 8 — > K is an entire solution to the minimal surface equation (MSE): 

V ■ ( , VF = 1 = in R 8 . (3.1) 



(3.2) 



The graph enjoys certain nice symmetries. Write x' — (u,v) S R 8 , where u, v £ R and u — 
v = \v\. Then F satisfies 

F is radially symmetric in both u,v, i.e. F — F(u,v) 
F(u,v) = —F(v,u), 

so that F vanishes on the Simons cone 

S = {u = v} = {x\ + ■ + x\ = x\ + ■ + xQ C M 8 . 

3.1. Geometry of the unit-width band around T. We will use powers of 

r(x) = \x'\ for x = (x',x 9 ) e M 8 x R = R 9 . 

to measure the decay rate of various quantities at infinity. As mentioned in the beginning, we are 
interested in the domain of definition for the coordinates (2.2): 

R 9 3 x ^ (y,z) eT xR x = y + zv(y), 

where v{y) is the unit normal to T at y € T such that v{y) ■ eg > 0. In particular, we are considering 
a type of domains which is a thin band around the graph T: 

B r (d) = {x e R 9 : dist(x,r) < d}. 

Lemma 3.1 (cf. Remark 8.1 in [ ]). There exists a small enough d > such that the coordinates 
(2.2) are well-defined in Br(d). 

The level surfaces for the signed distance to T 

T(z) = {x e R 9 : signed dist(x, T) = z) 

are prominent in our analysis as it will be necessary to estimate derivative operators on T(z), for 
small z. According to Lemma 3.1, the coordinates (2.2) are well-defined in a thin-enough band Br(d). 
Equivalently, the orthogonal projection onto T is well-defined in Br(d): 

7Tr : Br(d) — > T 7Tr(x) = y 
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and provides a diffeomorphism between T(z) and T for each z, \z\ < d. Thus, one can identify 
functions / defined on T with functions / defined on T(z) via the projection 7Tp: 

/ = /o7r r . (3.3) 

We will use the following lemma, quantifying the proximity of the gradient and the Laplace-Beltrami 
operator on T(z) acting on / to the gradient and the Laplace-Beltrami operator on T acting on /. 

Lemma 3.2 (cf. §3.1 in [12]). There exists < d < 2 small enough so that 

• The coordinates (2.2) are well-defined; 

• // we fix z with \z\ < d and assume f € C 2 {T{z)) and f € C 2 {T) are related via (3.3), we 
have the following comparison of their gradients, viewed as Euclidean vectors: 

|V rW /-(Vr/)o7^|=o(*^^)o7rr, (3.4) 
while the Laplace-Beltrami operators on T(z) and T are related by: 

W=(A r / + o(^ + ,J^))o Tr . (3,, 
Recall that we are interested in the decay rate of the quantities 

8 

where {fcj}| =1 are the principal curvatures of T - the eigenvalues of the second fundamental form. 
We prove the following useful estimates. 

Lemma 3.3. The k-th order intrinsic derivatives of the quantity Hi are bounded by 

\D h rMy)\ < 1 + % {y y (3-6) 

for some numerical constants Cm > 0. 

Finally, we would like to investigate how scaling space x — > a~ l x affects the estimates in Lemmas 
3.2 and 3.3. For a function / on T, define f a to be the corresponding function on T a = a _1 r: 

fa(y) = f(<xy)- 

Also, denote {kf }| =1 to be the principal curvatures of r Q and 

8 
«=1 

Lemma 3.4. Scaling space x — > oT x x has the following effects: 

• The intrinsic k-th order derivatives of f a , k = 0, 1, 2, . . . , scale like 

D$J a {y) = a k (D*f)(ay) y e T Q ; (3.7) 

• The quantities 

D h r H La {y) = a k+ \D$Hi)(ay) - O ( ' —-^-^j-t:) yeT a ; (3.8) 



1 + (ar(y)) l+k _ 

• If the coordinates (2.2) are well-defined in a bandBr(d), the coordinates 

R 9 3 x -t (y, z) G T a x R x = y + zv a (y) 1 

where v a {y) = v(ay) is the unit-normal to T a at y, will be well-defined in the band Br a {d / a) . 
Thus, the orthogonal projection onto T a is well-defined in Br a {d/a) : 

fr„ : Br a (d/a) -> T a irr a (x) = y; 



A FREE BOUNDARY PROBLEM INSPIRED BY A CONJECTURE OF DE GIORGI 



9 



If \z\ < d/a is small enough and f a £ C (T a (z)) and f a € C (T a ) are related via f a — 
fa 07r r a , the estimates corresponding to (3.4) and (3.5) take the form 

|Vr oW /a - (Vr./a) o ttpJ = ( aZ \ D f Ja \ ) on ra , (3.9) 

V 1 + ar / 

A ra(z) f a = ( A rJa +o(az l ^JA + ^ X ]^M\\ ^ (3J0) 



w \ V 1 + ar 1 + (ar)" 

Let us now turn to the proofs of the aforementioned lemmas. 

We will take advantage of the following local representation of the minimal graph T. At each 
y — (x' , F(x' )) € T denote by T = T(y) the tangent hyperplane to T at y. A simple consequence of 
the Implicit Function Theorem states that T can locally be viewed as a smooth (minimal) graph over 
a neighbourhood in the tangent hyperplane T. Concretely, if {ej}| =1 is an orthonormal basis for T 
and eg = v(y) is the unit normal to T, there exists a a small enough a = a(y) and a smooth function 
G :TC\B a {y) -> K so that in a neighbourhood of Xq, 

8 

{x',F{x')) = (x' ,F(x' )) + J2 t ie i + G{t)e 9 V\t\ < a. (3.11) 

i=i 

Moreover, G{t) satisfies the MSE H[G] — in {\t\ < a}. In [12] the authors establish the following 
key estimates for G, which provide the basis for the lemmas, stated above. 

Lemma 3.5 (cf. Proposition 3.1 and §8.1 in [12]). Fixy G T and let p = l + r(y). There exists a con- 
stant P > 0, independent of y, such that the local representation (3.11) is defined in a neighbourhood 
{\t\ < a(y)} C T with a{y) = j3p. Moreover, 

|AG(*)I<— i \D^G(t)\<^ m \t\<[3p (3.12) 
p p K 1 

for k € N and some numerical constants c,Ck > 0. Also, the unit normal v to T doesn't tilt signifi- 
cantly over the same neighbourhood: 

\u(t,G{t))-^y)\<— 1*1 <fip. (3.13) 
P 

The proof of Lemma 3.5 is based on Simon's estimate for the second fundamental form of minimal 
graphs that admit tangent cylinders at infinity, [20, Thm.4, p. 673], 



and employs standard MSE estimates applied on an appropriate rescale of G. 

Lemma 3.1, the possibility to define the coordinates (2.2) in a thin enough band around T, is a 
corollary of (3.13). 

Proof of Lemma 3.1. Assume the contrary: that there doesn't exist a d > for which the coordinates 
(2.2) are well-defined in Br(d). The coordinates will fail to represent a point x £ Br(d) uniquely 
when there exist two points i/i ^ j/ 2 £ T such that 

\x - yi | = \x-y 2 \ = dist(a;, T) < d, 

i.e. if 

x = yi + zv{yi) = y 2 + zv(y 2 ). 
We have \yi— 2/2 ] = I^H^Cj/a) - v {Vi)\ 2d. This means that if d is sufficiently small (d < | J 8(l+r(yi)), 
for example), y 2 lies in the portion of T which is a graph over T(yi) D {\t\ < (3(1 + r(yi))} - a 
neighbourhood of the tangent hyperplane at y\. But then (3.13) gives us 

cd 

\yi - 2/2! < d|^i - v 2 \ < — — 7 — ^- |t/i - 2/2 1 < cd\y x - y 2 \ 
1 + r(yi) 
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which is impossible whenever d is small enough so that cd < 1. □ 

Proof of Lemma 3.2. Let y 6 T(z), y — irr(y) and let T = T(y) be the tangent hyperplane to T at 
y. Use the Euclidean coordinates * on T to parameterize T(z) near y: 

i-> (t,G(*)) + 2i/(y(*)), 

where v{y{t)) = -^===a L is the unit normal to T at y(t) — (t, G(t)). Note that in these coordinates 

/(*) = /(*)• 

As before, set p = 1 + r(y) and use Einstein index notation. Because of (3.12), the metric tensor 
g{z) on r(z) computes to: 

gij(z) = Sij + GiGj + z(e-i + G,f) • djV + z(e.j + Gju) ■ div + z 2 diV ■ djV 
= gij (0) - 2zG lJ + z 2 d t v ■ d jV = 9ij (0) + 0{zp- 1 ) + 0{z 2 p- 2 ) 

while its inverse 

g»(z) = g i l(0) + O(zp- 1 ) 

for \z\ < d small enough. Noting that gij — <7y(0) is the metric tensor on T in ^coordinates, we see 
that the difference of gradients, viewed as vectors in Euclidean space: 

|Vr W /(»)-Vr/(tf)| = 
= \g»(z)d j f(0)(e i + <^(0)e 9 + zdns(Q)) - 3 «(0)9 J /(0)(e i + Gi(0)e 9 )\ = 
=0{z\D r f(y) Ip' 1 ). 
The derivatives of the metric tensor satisfy: 

dk[gij(z)] = dk9i 3 + 0(zp~ 2 ) 



so that 
Thus, 



dk[g i3 {z)] = -g U (z)d k [g lm (z))g^(z) = d^+Oizp- 1 ). 



A 



f(y) = -j^fdi (g ij C*) VWjldj /) 



A r /(») + O^p- 1 !^/^)!) + 0(^- 2 |£» r /(y)|). 



□ 



Proof of Lemma 3.3. Set p = 1 + r(y) and write the metric tensor of T around ?/ in the coordinates 
t, (3.11): 

Ay = 6 i:i + GiGj = 5 l3 + 0{P 2 ) |t| < /3p. 

Its inverse takes the form 

9 ij = Sij - = h + 0(/? 2 ) |*| < PP- 

Mind that constants in the O-notation are independent of p. Taking into account (3.12) we see that 
for |*| < Pp, k = 0,1,2... 

\D$ gij (t)\ < % |aV j '(*)I < 4- ( 3 - 14 ) 
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An easy consequence is the fact that the intrinsic k-th order derivative of a function / on T at y will 
be majorized by the i-derivatives of / at t = up to order k as follows: 

k 

i^/(»)i<E^ _fc i^/(°)i ( 3 - 15 ) 

for some numerical constants Cfc > 0. Since gij = 0(1) and g % i — 0(1), it suffices to show that 

Vj/ = dif + J2 c j(t)d.jf, where 

\J\<k (3.16) 
\D?cj(t)\ = 0(p- k+ \ J \~ m ) \t\<Pp, m = 0,l,... 

Here, of course, 7,7 denote multi-indices (e.g. if I = (ii,i%, ■ ■ ■ ,ik), |7| = k), V denotes covariant 
differentiation and 

V I f = (V k f)(d il ,d i2 ,---,d ih ). 

We'll prove (3.16) by induction on k = |7|. When k = 1, the statement is obviously true and assume 
it holds up to k — 1. For convenience, define the following transformation on multi-indices of length 
fc-1: 

of (jl,j2, ■ ■ • , jfe-l) = (jl, ■ ■ -,jl-l,j,jl+l, ■ ■ ■ ,jk-l) 1 < I < k - 1. 

So, if |7| = k and we write 7 = (ii, 7'), 7' = (12, ■ ■ ■ , ik), the covariant differentiation rule gives 

V I f = d il (Vrf)- J2 FU+^^f, 

l<Z<fc-l 
1<j<8 

where are the Christophel symbols for the metric tensor g in the coordinates t. We only need to 
check 7}™r^- = 0(p~ 1 ~ m ), which follows immediately from (3.14). The induction step is complete. 
Let us apply (3.15) to the second fundamental form 

(we use the Einstein index notation again). Observe that its i-derivatives decay like 

IA"W')| - o( P -^- m ) |t| < pp. 

Since 77, = Trace(LVf), 

D?H X = ■nacepn^'] 1 ) = 0(p- m ~ l ) \t\ < pp. 
Hence, (3.15) implies the desired 

|iW(tf)l < 



ifc V.t.W s- ^kl 
n l + k ' 



□ 



Proof of Lemma 3.4- The results in the lemma are obtained after simple length-scale considerations. 
Equation (3.7) is immediate. The principal curvatures kf(y) scale like distance^ 1 , so that 

K(y) = ak i( a y) ye r «. 

and thus 

T7 ; , Q (y) - a l H t (ay) = a'(77 ; ) Q (y). (3.17) 
We invoke (3.7) and (3.6) to obtain the full estimate (3.8). 
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We are left to check (3.9) and (3.10), which follow from (3.4) and (3.5) right after we note that 
(Vr Q ( z )/ Q )(y + zv a {y)) = a(Vr(az)f)(oc(y + zu(y)) 
(Ar a ( 2 )/o)(» + zv *{y)) = a 2 (Ar( tt2 )/)(a(j/ + ^(j/)), 
where / = / o n-p is the lift of / onto T(az). 

□ 

3.2. Proximity between V and the model graph. A more refined knowledge of the asymptotics 
of r is needed in order to carry out the construction of a supersolution to (1.1). To extract better 
information about geometry of T at infinity, Del Pino, Kowalczyk and Wei [12, §2] introduce a model 
graph Too, which has an explicit formula and which approximates T very well at infinity. Namely, 
the model graph 

r oo ={(x',x 9 )€M. 9 :x 9 =F 00 (x')} 
where : R 8 -> K solves the "homogenized" MSE: 



|VF„ 



0, (3.18) 



has the same growth (~ r 3 ) at infinity as F and shares the same symmetries (3.2). This determines 
the function F^ uniquely up to a multiplicative constant: if we use polar coordinates to write 

u = r cos v = r sin 9, 

the function takes the form Foo(r, 6) = r 3 g(6) where g{9) £ C 2 [0,7r/2] is the unique (up to a scalar 
multiple) solution to 

2W(20) U*n>W\' « 

,'(|) = o. 

which is odd with respect to 6 = 7r/4. For concreteness, we pick the g{9) which satisfies in addition 

<?'(!) = i. 

Del Pino, Kowalczyk and Wei then prove the following result quantifying the asymptotic proximity 
between the BdGG graph and the model graph. 

Theorem 3.1 (cf. Theorem 2 in [12]). There exists a function F — F(u,v), an entire solution to 
the minimal surface equation (3.1) which has the symmetries (3.2) and satisfies 

F OD <F<F 00 + — minji^, 1} m £[-,-], r > R (3.19) 
r 4 I 

for some constants C, Rq > and < er < 1 . 

The proximity of T and Too at infinity allows one to approximate geometric data and geometric 
operators defined on the non-explicit T with their counterparts on the explicit Too. To put it more 
concretely: one can use the orthogonal projection 7ir onto T to identify functions / defined on T with 
functions defined on far away from the origin: 

fco=fon r r m n {r > R} (3.20) 

for R > large enough. Then one can compare (Vr/) ° T to Vr^/oo and (Ar/) ° 7rr to Ar^/oo- 
Also, if we denote i}j =1 to be the principal curvatures of T^, and 



1 = 1 
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one can use the explicit H^ i to approximate the curvature quantities Hi associated with T. The 
ultimate goal is to approximate the Jacobi operator on T, 

Jrf = (A r + H 2 )f 
asymptotically with the Jacobi operator on Too, 

Jv^foo = (Ar^ +-Hoo,2)/oo- 

On the basis of Theorem 3.1, Del Pino, Kowalczyk and Wei establish the following list of results. 

Lemma 3.6 (cf. §8.2, 8.3 in [12]). Far enough away from the origin, r > R, for some constant 
< a < 1. 

• One can express T^ locally as a graph of a function Goo(t) over a neighbourhood of the 
tangent hyperplane T = T(y) to y £ T with r(y) > R. Moreover, for some constants Cj, > 0, 
k = 0,1,... 

\D k t {G-G QO )\\< Cl 



lt = 1 — r (yjk + l + <J ' 

where (t,G(t)) is the local parametrization (3.11) ofT. 

• The Laplace- Beltrami operator on T can be approximated with the Laplace- Beltrami operator 
on Too as follows: 

(A r J) o 7rp = Ar^/co + O (r"^^/^ + r^A^/col) , (3.21) 

where f and f^ are related via (3.20). 

• The quantities H 2 — \A\ 2 and H3 ofT are approximated by -ffoo^ an d Hqc.3, respectively, as 
follows: 

H 2 on r = H^ 2 + 0{r- A - a ) (3.22) 

H 3 on r = H 00 . 3 + 0(r- 5 -°). (3.23) 

• Therefore, the Jacobi operators on T and T^ are related by 

(Jr/)ovrr - Jv^foo + O (r^lD^U] + r-^l^/ocl + r^lf^) . (3.24) 

4. Construction of the super and subsolution. 
4.1. The ansatz. Define the L°° weighted norms 

\\f\\kx™ ( n ) = \\{l + r{y) k )f(y)\\ L ~ ( n ) 

for regions f2 C T. Use the short-hand || • \\k.oo when f2 = T. 

Recall that for a > small enough the coordinates (2.1) are well-defined in the band £>r Q = Br a (2). 
We will work with the following ansatz w : Br a — > R: 

w(y, z) = h%{y) + zK{y) + z 2 h«(y) + z 3 h^(y) + z 4 h«(y) + z 5 hUv), (4.1) 

where hf — hf (y) are functions on Br a , independent of the z- variable. The coefficients h",h 3 , h§ 
are explicitly specified in terms of geometric quantities associated with r Q : 

h? = l-\ A f + K K = "^(Ar Q + \A a f)\A a f ^ 

hi = \{\A a f + h?) hi? = \(A r JA a f - \AX) (4.2) 

n 5 — TTTT \ — n M4,a 7. 
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According to (3.8) of Lemma 3.4, the size of ft — (ft" — l),h% and their covariant derivatives up to 
second order on T, 

2 3 4 

h = 0[ - ,A \Dr a h\ = 0( - -. .A \D 2 h\ = o( - -. A , (4.3) 
Vl + (ar) 2 / 1 a 1 VI + (or) 3 / 1 la 1 Vl + (o;r) 4 / v y 



while for ft = ft'" ft'" ft" 



l + (ar) 2 /' 1 iQ 1 Vl + (ar) 3 /' 1 1q 1 \l + (ar) 4 , 

1 ) "3 i "-5 

" = °(iAA '^-^(AA I ^'" = (TTW)' ( "> 

We set ft" = 0. The coefficients ftp and ft 2 will be specified later so that the ansatz meets the 
supersolution conditions in Definition 2.2, but from the very start we will require that they satisfy 
the following properties: 

• ftg > is strictly positive and scales like 

h%(y) = aPh (ay), 
where < p < 1 and fto £ C 2 (T) is positive with 

H^rMs.oc + ll^rftolkoo + HMi.oo < C x . (4.5) 
for some positive constant C\. So, 

/ nP \ / rr 1+p \ / rv 2+p 

K = o( T ^ 7 -A, \Dv a K\ = o(—— — ), \D 2 K\ = o( 



.l + (or)/' 1 " Ul Vl + (ar) 2 /' 1 1 " ul Vl + (ar) 3 
• ft 2 equals 

h%=\{\A a \ 2 h% + hf), 
where the correction ft 2 " scales like 

h' 2 a (y)=a 2+ Ph' 2 (ay) 

and ft 2 £ C 2 (r) is positive with 

IPr^lkoo + ll^rft 2 ||4,oo + ||ft 2 ||3,oo < C 2 , (4.6) 
for some positive constant C 2 . Thus, 

K = 0[- ^A \D Va h%\=0[- --A \D 2 h 2 l \=0[- -—,). 

1 Vl + (w) 3 / 1 2i Vl + (or) 4 / 1 La 21 Vl + (ar) 5 / 

Remark 4.1. At first look, the choices for hf above may seem somewhat arbitrary but they are 
prompted by the supersolution conditions. The fact that we expect the solution to behave asymptotically 
like z suggests that hf m 1 to main order. Thus, the main order term in H ra ( z )d z w is z\A a \ 2 which 

has to be cancelled by the z 1 -term in d 2 w: thus, h% rs . Now d z w ~ ft" + z 2 A^ and since w 
achieves values ±1 at z ~ ±1 and |Vlo| ss d z w, the supersolution gradient condition demands that 
we refine ft" to equal ft" w 1 — A^ . The form of h 2 is contingent upon the fact that w(y, •) attains 
the values ±1 asymptotically at z± ~ ±1 — ftp , so that 

d z w(y,z ± )^l±(2h 2 <~\A a \ 2 h%), 

requiring the positivity of ft' 2 " = 2ft£ - \A a \ 2 h%. The remaining choices ( and further refinements ) 
are made so that no terms that decay at a rate r~ 4 (and no better) at infinity are present in the 
expansions of Aw or d z w(y, z±). 

All this will become transparent once we carry out the computations of the Laplacian of w in 
Lemma 4.1 and of the gradient of w on {w — ±1} in Lemmas 4-2 and 4-3 below. 



A FREE BOUNDARY PROBLEM INSPIRED BY A CONJECTURE OF DE GIORGI 



15 



NB. In what follows the constants in the O-notation depend solely on p, Ci, C2 and the minimal 
graph r, but not on the scaling parameter a. 

Lemma 4.1. The Laplacian of w in Br a can be estimated by 

Aw(y, z) = (A Va{z) + \A a \ 2 )h° + hf - z 2 H 3 . a + o( — — j . (4.7) 
Proof. Compute in succession: 

d z w = 1 - JM + h >° + z (\A a \ 2 h« + hf) + z ^ 2 + h " + 5z 4 ft« (4.8) 



9> = \A a \ 2 h% + ti 2 a + z(\A a \ 2 + hf) + 20z 3 /£ (4.9) 

v5 



H Va{z) d z w= (z\A a \ 2 + z 2 H^ a + z 3 H^ a +0(z 4 - + ° ar)5 )) d *™ 



, - lAf + o( TT ^)) + 2 = ( ff 3,„ + o( T |^ p ) ] («o> 
^(^ + *, o + 0( r " 6+ ' 



2 4 



+ (ar) 6 // V 1 + (or) 5 
Because of (3.10) and (4.4), 

\Ar aM hT\+\^r aM hT\ + |A r „ w fc?| = o( T 

Ar aM ^ = 



+ (ar) £ 



1 + (ar) 5 
Ar Q(i) |A Q | 2 = A r JA Q | 2 + 0(z r 



+ (ar) 1 



so that 



Ar. W - - - ^ + **f£ + O(^). (4.1, 



Combining (4.9), (4.10) and (4.11) we derive that in Br a 

Aw = (A Fa(z) + \A a \ 2 )h% + h' 2 a - z 2 H^ a 

+ z 3 (^f^ + 2, K -^- Hi 



(A Ta(z) + \A a \ z )h% + h' 2 a - z l B^ a + 0[j 



+ (ar) 5 J 



a 



Now we would like to determine how far the level surfaces {w = ±1} stand from the graph T a . 
Note that for \z\ < 2 and uniformly in y e T a , we have w = z + 0(a p ) and d z w = 1 + 0(a 2 ). Thus 
for all small enough a > 0, w(y, ■) is strictly increasing and attains the values ±1 for unique z±(y) 
with \z ± (y)\ < 2. 

Lemma 4.2. For all small enough a > (so that z± is well-defined), 
d z w(y,z ± (y)) = l±h' 2 a - 1 ° j -ti?K + 0( 



1 + (ar) 5 
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Proof. Fix y £T a and let us estimate z±{y). Write z± = ±1 + 6±, where S± = o(l) as a — > 0. We 
compute 

w(y, z ± ) = ft? + (1 - ^A(±l + 5 ± ) + ft?(±l + <5±) 2 + 



2 



14 



2 



-(±i + s ± y + o 



Thus, 



6 v ; Vl + (ar) 4 

= ft«±(l-JM) + ft-±JM + 

+ ( y ± (i±2/ l2 a ) + 4(±^ + ^) + 4^£ + O u + w 



+ (ar) 4 

so that 



5± = o(ft? 



ft? + = O 



which in turn implies 



\A I 2 / a 2+2p \ 

fe— «± L f L -** +0 (iw)- (412) 



We can now estimate d z w(y,z±(y)): 



d zW {y,z±(y)) = i - ^ + fti Q + 2ft? (±i - ft?)+ 

\A a \ 2 + hf \A a \ 2 2 a 

+ 2 (±1 "o± 3 J + 5 ^5 +°Vl + («r) 

= 1 + (±2h% T |A| 2 ft?) + (-2ft? ft? + 

h' a I A I 4 / ™ 2 + 3 P 

+ (ft'i Q + % + ^+5ft?) + 



. 1 + (ar) 5 
j_ ± h ,« _ IA| 2 (ft?) 2 _ h ,a h a + q 



1 + (ar) 5 



v 4 +P 



□ 



Straightforward derivative estimates using (3.9) yield 
Lemma 4.3. We have 

|v rQ(2±)W | 2 = |v ra(2±) ^| 2 + o( TTW ) 

and thus 

|Vu;| 2 (y, z±(v)) = (d z wf + |Vr Q(2±)W | 2 = 

1 + (ar)5 j 

Our ansatz has the very nice, extra feature that it is strictly increasing in Br a in the direction of 

eg. 

Lemma 4.4. For all a small enough 

d Xg w > in Br a ■ 



A FREE BOUNDARY PROBLEM INSPIRED BY A CONJECTURE OF DE GIORGI 



17 



Proof. Computing in the coordinates (2.1) 

d X9 w(y, z) = V«i • e 9 = (V ra ( z) w + {d z w)v(y)) ■ e 9 

> — / vr (z)UJ\ 



y/l + \VF a \* Vl + (ar) 2 

Since 

1 

> 



1 + OM (t ^). ,4,3, 



v/i + |v^ Q | 2 " i + K) 2 

for some positive constant c > (see Remark 8.2 in [12]), (4.13) yields 

d xg w > 

for all small enough a > 0. □ 

4.2. Supersolutions for the Jacobi operator. As we have noticed from Lemma 4.1, the sign of 
Aw depends crucially on whether the Jacobi operator 

Jr Q = A r „ + \A a \ 2 

admits positive supersolutions that satisfy appropriate differential inequalities. We will show that the 
Jacobi operator Jr on the (non-rescaled) minimal graph L admits the following two types of smooth 
supersolutions: 

• Type 1 is a positive supersolution h € C* 2 (r) such that for some < e < 1 

Jrh(y) < ' 



1 + r 4 + e (y)' 

• Type 2 is a positive supersolution h G C 2 (r) such that 

x hU , w l%)-^/4| 

1 +r 6 {y) 

The Type 1 supersolution is readily provided by [ , Proposition 4.2(b)] (our Proposition 4.1 below 
is a straightforward modification). We construct the Type 2 supersolution in Proposition 4.2 and the 
supporting Lemma 4.6. 

Proposition 4.1. Let < e < 1. There exists a positive function h E C 2 (r) such that 

\\Dj.h\\ 4+e ^ + ||I>r/i||3+ 6 ,oo + \\h\\ 

and 

Jrh < 1 7 — . 

1 + r 4 + e 

Proposition 4.2. There exists a non-negative function h £ C 2 (r) such that 

ll-Dr^lkco + \\Drh\\ 2 ,oo + \\h\\i,oo < 00 

and 

1 + r J 

Moreover, there is a ^ < r < | (e.g. t = ^) such that for every < 6 < 8' < | 

||'l||l+«r,£ 00 (S(-5')) + H-Cr^||2+(5r,L~(S(-(5')) + II A^ll 3+<5t,L°= (S(-<5')) < 00 ( 4 -14) 

w/iere 5(-<J') = {|0 - f | < (1 + r)- s '} C T. 
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Before we venture into proving these two propositions, recall that 

Jrh < f in the weak sense if 
(Jrh - f)[4>] < for all non-negative € C*(r). 

Above we have used the notation 

fl<t>] = ! f<f> for / g ii oc (r) 



(J r /i,)[0] = j -\7 r h-\7 r cf>+\A\ 2 



where test functions 4> € C*(r). 

Let us make the important remark that the operator Jr satisfies the maximum principle. 

Remark 4.2 (Maximum principle for Jr)- Since ho := ^ 1+ | VF p - > solves Jr/io = (see (A. 2),), 
the elliptic operator 

L := h A r + 2V r /i • V r 

satisfies 

J r h = L(h/h ). 

Thus, if h is a supersolution for Jr (in the weak sense) in a bounded domain U C T and h € C(U), 

> Jrh = L(h/h ) in U, 

so that the quotient h/ho doesn't achieve its minumum at an interior point of U unless h/ho is 
constant in U. 

In fact, we will construct the supersolutions in Propositions 4.1 and 4.2 as solutions to appropriate 
elliptic differential equations rather than inequalities. This approach will pay off, because in the end 
we will automatically possess global smooth supersolutions, whose first and second derivatives will 
have the appropriate decay rates at infinity. Specifically, we will investigate the linear problem 

J r h = f in T, (4.15) 

where h and / are in appropriately weighted Holder-type spaces. As usual, we first study the problem 
(4.15) in bounded domains Tr :=m{r< R} 

Jrhn* = f in T R 
h B = on dT R . 

Because Jr satisfies the maximum principle, the problem (4.16) is uniquely solvable for all R. In 
order then to run a compactness argument which takes a sequence /i_r„, R n /* oo and produces a 
globally-defined h : T — > R that solves (4.15), we need two important ingredients - the existence of 
suitable global barrier functions and an a priori estimate (Lemma 4.7) for the solution to (4.16). 

We first exhibit functions that are (weak) supersolutions for Jr far away from the origin. Later, 
we will be able to modify and extend them to barrier functions on the whole of T. 

Lemma 4.5 (cf. Lemma 7.2 in [12]). Let < e < 1. There exists a positive function h, such that 
for some R > and constants c, C > 

JvKy) < -Y^A+i m Mv)>R} (4.17) 

7^72+7 - h ( y ^ - T+72T? m {r(v)>R}- (4-18) 



A FREE BOUNDARY PROBLEM INSPIRED BY A CONJECTURE OF DE GIORGI 



19 



Proof. It follows from the existence of the Type 1-supersolution hi j0 o G C 2 (roo) for Jr^ far away 
from the origin {r > R} (See (A. 7) of Appendix A): 

Jr x h 1 ^(y)<- i + r 1 ^ )4+e y€T 00 n{r>R}. 

Use the orthogonal projection ir-p to lift /ii l00 to a function hi on L n {r > R} 

hi o 7rr = /ii.oo on r n {r > i?}. 
Then according to (3.24) and the gradient and hessian estimates in Lemma B.l (see Appendix B), 
J r h(irr(y)) = Jr ao hi t00 (y) + 

+ O (r- 2 -°\D^h hoo \ + r- 3 -°\D Too h hoo \ + r~ 4 ~°\h hoo \) (y) 
1 

< -- 



2(l + H+«(jfl) 

for r(y) > R large enough. Equations (4.17) and (4.18) are obtained once we note that, according to 
Lemma 3.6, y := 7Tr(y) is very close to y : 

\y-y\ = 0{r- 1 -°{ y )) 

in {r > R} for a large enough R. 

□ 

Lemma 4.6. There exists a locally Lipschitz, non-negative function h, which is a weak super 'solution 
for Jp away from the origin and which satisfies 

Jrh(y)<- 16 ~^ 1 on {r(y) > r } (4.19) 
1 + r A 

for some large enough r > 0. Moreover, 

h = (\ e -*W + l) (4.20) 
V l + r l + r 2 + e J V ; 

for some t£ (§> g) an d some e G (0, 1) (e.g. r = | and e = | do the job). 

Proof. The construction of the weak supersolution in this case is achieved by patching up two smooth 
supersolutions, defined on overlapping regions of T, via the min operation. The resultant function is 
obviously locally Lipschitz. 

One of the building blocks is the Type 2 supersolution for Jp^ at infinity (A. 8)- call it /i ext G 
C^Too n {tt/4 < 6 < tt/2}) here: 

where q2(0) has the following expansion near = ?: 

= (A - J)>o + a 2 (^ - ^) 2 + • • • ) a > 0, 

and 

TT73~~ 0e( 4'2- 

for some large r > 0. Define for 

- 2 < a 2 < ax < (4.21) 
the following subregions of the model graph L^ 

Too.int - {\B - ~| < n {r > r } G 
Too.oxt = {\6 - \\ > r a ?} n {r > r } C r M . 



Jr^feext < — , 3 0e( T ,-) and r > r . 
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Note that L^nt and r^^t have a non-empty overlap and that they cover all of n {r > r^}. 
Define km G C 2 ^) by 

- (9 - Tr/4) 2 r 3+s 



Vl + IVFool 5 



(4.22) 



for some < S < 1, which will be specified later, and extend /i ex t on the whole of f]{r > tq} so 



that it is even about 9 = ?: 



7T 

h c *t(r,0) = h cx t{r, - 



The goal is that the lifts of h[ n t, h cx t onto L: 

^int(7rr(y)) = ^int(jf) h C Kt(^r{y)) = hext(y) ye 

corrected by an appropriate asymptotic supersolution of Type 1 (given by the previous Lemma 4.5), 
will satisfy the desired differential inequality (4.19) in the respective regions 

Tint = Trr(rooant) C T and r oxt = 7rr(roo,ex±) C T. 

Applying the gradient and hessian estimates of Lemma B.l from Appendix B and the fact that 
—2 < a2 < a%, we derive that in r m -int 

\knt\ = 0(r 1+s + 2 ^) 

\D rao knt\ =0{{9- ^) 2 r s +r- 2+s \6- ~|) = 0(r s+2a ' + r^ +s - 2 ) = 0(r s+2ai ) 
iD^h^l = o((0 - ^) 2 r- 1+S + r- 3+s \6 - J + r- 5+i ) = 0(r- 1+5+2Ql )■ 
On the other hand, h cxt satisfies in r^ext 
\h**\ = O(\0-~\ r r- 1 ) 
\D r Jh cxt \ =0{\9- \Vr- 2 + r" 4 |0 J- 1 ) = 0(\9 ^r~ 2 ) 
\D 2 r J cxt \ = 0(\9-\Yr-* + r-*\9-\r^r-i\9-\Y- 2 ) = 0(\9 - \Yv~ 3 ). 
Thus, the proximity (3.24) between Jp and Jp^ implies 
Jrh int (My)) = 0(r- 1+5+2a i(y)) y E 
Jrh e *(Mv)) < - m f + ~^ T + 0(\9(y) - ^/^r^iy)) 

l jm-jr/W C|fl(y) - tt/4| 
S 2 l+r 3 (y) S l + r 3 (y) 2^ W*t 

for large enough r(y) > r^. According to Lemma 3.6, if y — irr(y), 

l»-»l = 0(r(»r 1 " ff ) 

for some < a < 1 and r(y) > r$ large enough. Therefore, the pair (r(y),9(y)) is asymptotically 
equal to (r(y),9(y)): 

\r(y) - r(y)\ = 0(r- l -°(y)) \9{y) - 9(y)\ = 0( r - 2 ^(y)). 

Thus, 

- IT Ml I0(V>-"7tMI 

0(r(y)- 



_ i%)-V4| , ,„ , , Ti 

l+r 3 (y) l + r 3 (y) 
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so that 



= 0(r- 1+s+2 ^(y)) yeT iat (4.23) 



Jrh' <- — r>r Q (4.25) 



Let hi be the supersolution for Jr, provided by Lemma 4.5: 

1 

i +7 4 

for some < e < 1 which we'll pick shortly. Below we will define the functions h\ and hi, which will 
be supersolutions for Jr on Tint and r ext , respectively, and patch them into a (weak) supersolution 
h, defined on T n {r > r }, via the min-operation: 

h = min(/ii, hi). 

In order for the operation to succeed, we have to verify the following: 

• hi := h mt + h! satisfies the differential inequality (4.19) in Tint for large ro- 

Mh iat + h') + ~ ]9 ~ n/ ^ <0 in r int . 
Because of (4.23) and (4.25), it suffices 

3 6 S 

-1 + S + 2ai < -4 - e a\ < 

• hi := /lext + ^' satisfies (4.19) in r ext for large tq: 



Jr(/W + h!) + - 1 1 | 7 3 1 < in y e r cxt . (4.26) 



C ]6>-tt/4 ] 
2" 1 + r 3 

By (4.24) and (4.25) this holds for a sufficiently large ro- 
• hi < hi on Tint \ r cxt and hi > hi in r cxt \ Tint, i.e. we would like to have h m t < ^cxt on 
Tocvnt an d ^int > ^ext i n ext . This will be the case for large enough r if 

2 -\- 5 

a 2 < < ai (4.27) 

Z — T 

Collect conditions (4.21), (4.26) and (4.27) in 

2 + ^ -3-e-S ,„„„. 
-2<a 2 <-- < ai < . (4.28) 

Z — T Z 

Moreover, (4.28) needs to be compatible with 

, , , 1 2 
e,<5e(0,l), 3 <T< 3- 

Condition (4.28) is fairly tight, but not void: indeed, for 8 = |, e = | <E (0, 1) and t = | 6 (|, |), 
we have 

20 29 

-2 < a 2 < < oti < . 

~ 11 16 

So setting the parameters appropriately, we can conclude that 

{hi in r int \ r oxt 

mm(/ji,/* 2 ) in r int n r cxt (4.29) 
hi in r oxt \ r int 

is a weak, locally Lipschitz, supersolution for Jr in r > tq for a large enough ro > that satisfies 
(4.19) and (4.34). □ 
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The second ingredient is an a priori estimate for the solution h to (4.16). Introduce the Holder- type 
norms: 

/(l/l)"/(«2) 



\f\c(n) = sup 



dist r (yi,2/ 2 ) 7 

||/llk,CT(n) = ||/||fe,i~(n) + 11(1 + r k+ ' t )\f\c-i(c lH1+r) ( y )nn)\\L'^(n) 
where Q C T, k > 0, < 7 < 1, distr(yi, 2/2) is the intrinsic distance on Y and 

8 

Cr(«) = {^2uh + lv(y) :\t\<r, le R} 

i=l 

is the infinite right cylinder with a base B' r := {\t\ < r} C T y on the tangent plane T y to y G T (see 
(3.11) to recall notation). 

We will now establish the following regularity estimate. 

Lemma 4.7 (compare to Lemma 7.5 in [12]). Let R > be finite or infinite and assume h G C 2,1 (Y r) 
is a solution to (4.16) with f £E C^(T R ). Then 

Pr/IU+2,C7 ( r R/2 ) + ||£>r/IU+i,£~(r fl/2 ) < C(\\h\\k^(r R ) + ||/|| fe +2,c, ( r«)) (4.30) 
with a constant C > 0, independent of R. 

Proof. The proof is based on a rescaling technique. We may assume 

II^IU,i~(r H ) + ll/IU+2,C-r(r Ji ) < 1- 

Pick y € set p = 1 + r{y) and express the operator Jr in Cp p {y) PI r# using the coordinates t 

(3.11): 

= g»d%h + d^d 3 h + ^-^d,h + \A\ 2 h = 
2 \9\ 

= a^dfjhit) + VdMt) + \A\ 2 h(t) = f(t) t e B' Pp . 

Rescaling to size one, 

h(t)=p h h(py), f{t) = p k+2 g{pt), a^(t)=a ij (pt), V = pb\pt), 

we get 

a^dffhit) + b l d~h{t) + \A P \ 2 (t)h(t) = f{t) in B' p . 

Recall the standard Holder norm of a function q defined on a domain [/CR 8 : 

., ., „ „ \q(t) — q(s)\ 

\\q\\o(U) ■= h\\ l°°(u) + sup — — . 

t^seu \t-sp 

Because of the estimates (3.14) on the metric tensor g and its derivatives, and the estimates (3.6) on 
the second fundamental form \A\ 2 and its derivatives, we can bound 

ll« 4J llc,(Bp,||6 l ||c^),|||v4 P | 2 ||c,^) <^ 
by a universal constant K. Thus, by interior Schauder estimates, 

WDlhWc^B'^) + \\Dth\\Loo { B' m ) < C(IWU~(^ /2 ) + ll/1c(^ /2 )) < C 

so that 

p k+2+ -'\D 2 r h(y)\ CHrRnC ^ My)) + P k+2 \Dlh{y)\ + p k+1 \D r h(y)\ < C", (4.31) 
for each y £ T R / 2 - □ 
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Let us now show that 

J r h = f in F, 

where the right-hand side ||/|U+2.c^(r) < 00 f° r some k > 2, is uniquely solvable when ||/i||fc )00 < 00. 

Proposition 4.3. Let k > 2, 0<7<1 and \\f\\k+2,ci(T) < 00 • There exists a unique solution 
h G C 2 (r) to (4.15) such that ||/i||fe,i<»(r) < 00 • Moreover, 

\\D 2 T h\\ k+2 ^ (T) + \\D r h\\ fe+i,L°°(r) + ||^||fc,L°°(r) < C||/llfc+2,CT(r)- (4.32) 

Proof. Uniqueness follows from the maximum principle (Remark 4.2) and the fact that \h/ho\ < 
Cr 2 ~ k — > as r — > 00. 

To establish existence, consider the Dirichlet problem in expanding bounded domains: 

JtK = f in T Rn 
h n = on dT Rn , 

where R n 00. First claim that 

\\hnh,L~(r Rn ) < C\\f\\ k+2iC , { r Rn ). (4.33) 

for some constant independent of n. Assume not; then there is a subsequence (call it R n again) such 
that 

\\h n \\k,L°°(r Rn ) > ^||/||fe+2,CT(r Jin )- 
If we set /„ = f/\\h n \\ k ,L°°(r Rn ), h n = h n /\\h n \\k,L°°(r Rn ), we see that 

Jrh n = fn 

with ||ft-n|U,L° c (r ii?1 ) = 1 an d ||/n|U+2 ) CT(r H ) < \/n. The a priori estimate (4.30) implies, after 
possibly passing to a subsequence, that h n converge uniformly on compact sets to a C 2 (r)-function 
h with ||ft||fc )00 < co which solves 

J r h = in T. 

Uniqueness requires that h = 0. Let h'^ be the supersolution for Jp provided by Lemma 4.5 with 
some < e < k — 2 and ro large enough: 

1 c 

■/rfeoo <- 1 + r4+e i fc oo > x + r2+e r>r 0. 

Since /i„ — > uniformly on compact sets, s n := sup r ^ h n — > 0. Therefore, 

±/i„ + Unh'oc > on r = r and r = R n 
Jr{±h n + Unh'oo) < in r < r < R n 

for [i n = max{ SnC" 1 , ^} — > 0. An application of the maximum principle yields 

|fcJ<rtXo in r <r<i?„. 
Combine this with the fact that ||/in|U,£°°(r ro ) < s n r o *° conclude 

\\h n \\k,L°°(r Rn ) -> as n^oo 
which is a contradiction. Hence, (4.33) holds and the a priori estimate (4.30) becomes 

ll-C'r /l nlU+2,c^(r H?i/2 ) + llA^||k+i,i«>(rj W 2) + \\ h n\\k,L°°(r Rn/2 ) < C\\f\\k+2,c-r(r R „) 

for some constant C, independent of R n . Now a standard compactness argument produces a C 2 (r)- 
function /i which solves (4.15) and satisfies the estimate (4.32). □ 

Proposition 4.1 is an immediate corollary. 
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Proof of Proposition 4-1- Let h solve (4.15) with a right-hand side / = — 1+ * 4+ c ■ We are only left 
with checking that h is strictly positive. This is a consequence of the strong maximum principle 
(Remark 4.2) and the fact that \h/h (y)\ = 0(r~ e (y)) -t 0, as r(y) -> oo. □ 

Now we would like to construct a global barrier function (not-necessarily smooth) for (4.15) with 
a right-hand-side 

|0-tt/4| 



f = - 1 



1 



Lemma 4.8. There exists a globally defined, locally Lipschitz function h > which is a weak super- 
solution for Jr and which satisfies 

10 -7T/4I 

J T h<-- V- 1 in T. 

1 + r 3 



+ (4.34) 



Moreover, 

h = nfl 

l + r 1 + r 2+£ 
for some t £ (|, |) and some e € (0, 1) /e.g. r = | and e = |J. 

Proof. Let /i^ be the weak supersolution for Jr in r£ o , provided by Lemma 4.6: 

(JrC-/)M<0 

for every non-negative € Cg(r£ o ). Now let ^ g C°°(r) be a non- negative cutoff function such that 

tp(y) = for r(y) < ro and VKz/) = 1 f° r r (y) > ?*o + L 
Define a function /i" on the whole of T by 



fc"(y) 

Finally set 



in r < r 

^(y)h'L(y) m r > r o 
h = cti + h", 



where h' is the supersolution provided by Proposition 4.1 and C > is some large constant, to be 
fixed shortly. Now for any nonnegative 4> € C*(r), the fact that (J-ph 1 ^ — f)[ip<f>] < implies 

{Jrh - f)[<p] = C(J r ti)[4>] - f[4>] + J -^VrV ■ V r - V r C • (^V r 0) + |A 2 |^V0 
= C(J r ti)[cf>] - (1 - + (JrO^] - W<t>] 

+ [ (2V r ^ • Vrt& + C^r^ < -C(J r h')[<t>] + k[ct>] 



where k is a bounded function, compactly supported in r ro+ i. We were able to carry out the 
integration by parts, since h'^ is locally Lipschitz. Taking C > large enough we conclude that 
Jrh < / globally, in the weak sense. □ 

We now possess all the means to prove Proposition 4.2. 

Proof of Proposition 4-2. Pick / g C° n (T) such that / < and 

/o,r® = -M yeT oo n{r>r 0} 
l+r 6 (y) 

for a large enough ro. It is not hard to verify that ||/||3 ) cr(r) < oo by transferring the computation 
onto Too via (3.4) and employing the gradient estimate in Lemma B.l. 
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Let h n solve the Dirichlet problem (4.16) in the expanding bounded domains T Rn , R n oo 

Jrh n = f in T Rn 
h n = on dT Rn . 

Since / is non-positive, the weak maximum principle implies h n > 0. Let h! be a Type 1 supersolution, 
provided by Proposition 4.1, and let h" be the weak Type 2 supersolution which we constructed in 
Lemma 4.8. Noting again that 

\9(y)-w/4\ |%)-tt/4| 



l + r 3 (y) ' l + r 3 (y) 
for some a > 0, we obtain 



+ 0(r(y)- & - CT ) y = n r (y) 



J r (-h n + h" + Cti) < 

for a large enough C . Moreover, since —h n + Ch' + h" > on dT Rn , the maximum principle implies 

< h„ < h" + Ch 1 in T Rn . 

Thus, || /in || 1,1, ""(rR ) ^ C f° r an absolute constant C independent of n. We can now employ the a 
priori estimate (4.30) into a standard compactness argument that yields a non-negative C 2 -function 
h solving 

J^h = f in L 

with 

0<h< h" + O(h') 

,- - - (4.35) 

ll-Dr^lkoo + ||-Dr7i||2,oo + \\h\\i,oo < oo. 

After possibly correcting h by a supersolution of Type 1, h = h + ch' , we can conclude that 

Jrh(y) = f + cJrh'< J l V l-;f . 

l + r 6 {y) 

To establish the second statement in the proposition, namely the refinement of decay of h near 

3 



9 = 7r/4, we notice that on S(—6) with < S < 3 



\ 1 + r l + r 2+ / 

as St + 1< 2. Also, \\f\\ 3 +8T,Ch(S(-S)) < °°- 

Then an argument, based on rescaling and interior elliptic estimates - absolutely analogous to the 
one for the a priori estimate (Lemma 4.7) - gives us the interior estimate (4.14) (for h and thus for 
h itself) on S(-S') <s S(—8). There is a caveat: the same argument will carry through to the present 
situation, once we ascertain that S(—6) contains "balls" of size ~ r, centered on points in S(—6') far 
away from the origin. More precisely, we want for some r large enough, 

C My) (y) n L C S(-S) for every y e S{-6') n {r(y) > r Q } (4.36) 

Note that according to Lemma 3.5, the fact that T is a graph {(t, G(t)} over B'p r ( y )(y) with 

\G(t)\<Cr(y) 

implies 

Cf3r(y){y) n r c B cor{y) ( y ) n r 

for a large enough numerical constant c$ > 0. Suppose that (4.36) is not true: then there exist 
y' G dS(-S') and y 6 dS(—5) with r' = r(y'), r = r(y) arbitrarily large such that \y — y'\ < c r' . 
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Denote the projections of y' and y onto M 8 by (u',v') and (u,v), respectively. Obviously, r/r' ~ 1 
and for < S < 5' < 2 

|y' - y| 2 = |u - u'\ 2 + \v- v'\ 2 + \F{r', (1 + r'y 8 ' ') - F(r, (1 + r)- 5 )\ 2 
>(u- u'f + (v- v' f + \F(r', (1 + r'y 5 ') - F(r, (1 + r)- s )\ 2 
> (r' - r) 2 + iF^r', (1 + r')- g ') - F^r, (1 + r)- s )\ 2 - C'r^ 

r \ 3-5 



> -C'r- 2a + (/ - rf + c' 



2 (/) 2(3 ~ 5,) » (r 1 ) 2 



which is a contradiction. □ 

4.3. The free boundary super and subsolution. In correspondence with the form of the super- 
solution ansatz (4.1), define the subsolution ansatz v : Br a — > K by 

v(y, z) = -h«(y) + zh°(y) + z 2 (~h«(y)) + z 3 h«(y) + z 5 h%(y). (4.37) 

Since we require /iq > and h% > 0, we will automatically have v < w in Br a ■ Also, 

/ a p 

0<w-v = 2{K + z 2 h") = O 

v 11 \ 1 + ar 

Proposition 4.4. Fix < p < 1. There exist C\,C% > and ao > sttc/i that for all small enough 
a < ao, w given by (4.1) satisfies 

Aw < in Br a 

|Vw| 2 > 1 on {w = 1} and (4.38) 
|Vw| 2 < 1 on{w = -l], 

while 

Av > m 

|Vw| 2 < 1 on {v = 1} and (4.39) 
|Vw| 2 > 1 on{v = -l}. 

Moreover, < w — v < \, d Xg v > and d Xg w > in Br a ■ 

We immediately derive as a corollary: 

Corollary 4.1. Let i> , w, < a < ao be as in Proposition 4-4 above. Then the function W : M 9 — > M. 
given by 

!w(x) for x G B Fa n {\w\ < 1} 
1 for x G {B Ta H < 1}) C n {x 9 > F(x')} 
-l /or x g (Sr a n {H < 1}) C n {x 9 < F(x')} 

is a classical strict supersolution to (1.1), while the function V : M 9 — > R, given by 

( v(x) forx£ B Fa H {|u| < 1} 
F(x)=i 1 /or a; G (B Ta f~) {\v\ < 1}) C n {x 9 > F(x')} 
{ -1 for x G (S Fq n{|u| < l}) c n{x 9 <F(z')} 

is a classical strict subsolution. Moreover, < W — V < both V and W are monotonically 
increasing in xg and strictly increasing in xg inside f2j n (V), £li n (W), respectively. 
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Proof of Proposition 4-4- Fix some < S < \ and < e < St, where \ < r < | is provided 
by Proposition 4.2 and let h! > 0, h" > be the Jr-supersolutions given by Lemma 4.1 and 4.2, 
respectively. Remember that we only need to set the values of ho and h' 2 in order to determine the 
ansatz (4.1) completely. So, let 

ho = h + h , 

and 

, _ It 1 ccos 2 (26>)n 
2_ 2Vl + r 4 + e+ 1 + r 3 J' 
where c> is such that 2ccos 2 (26>) < \6 - 7r/4|. Set 

Cl = ||-D r ft,o||3,oo + ||£>r^o||2,oo + ||/lo||l,oo 
C'2 = HlMkoo + ll-Pr^lkoo + \\h' 2 \\ 3 ,oo- 

Claim that for all small enough a > 0, Aw < in Br a ■ This is a consequence of the following 
computations. 

• For a > small enough, 

tt/4| /a 3 |(9-7r/4| 



«/r — Hi n, < — — — — — — -|- ( 

" ' l + (ar) 4 + £ l + (ar) 3 V 1 + (arf 



a 2+ P a 



2 +p |0-tt/4| 



l + (ar) 4 + e 2 l + (ar) 3 ' 

so that 

Jv a K - Z *H 3 a + h>°<-^ ( + i^T) . (4.40) 
° u ' 2 ~ 4 Vi + C" 7 ") l + (or) 3 / v ; 

In S a (-1) = {\6 - f | < (1 + ar)- 1 }, Proposition 4.2 and (3.10) imply 

(Ap q(2) + \A a \^ = J Ta K + 0( 1 + K)4+ , r ) • (4.41) 

Then (4.7), (4.40) and (4.41) yield the desired 

Aw(y, z) < for y € and (y, z) € #r Q , 
and all small enough a > 0. 

In SS(-l) = {|0- f| > (1 + ar)" 1 }, (4.40) can be estimated further by 



Because of (3.10) we have 

(A ra(z) + |A Q | 2 )^ = Jr a hS + o( 1+ ). (4.43) 

Thus, (4.7), (4.42) and (4.43) yield 

Aw(y, z) < for y G S£ and (y, z) e Sr , 
and all small enough a > 0. 
To verify that necessary gradient conditions (4.39) are also met, we need to check that for small 
enough a > 0, h' 2 a majorizes both |A a | 2 (/ig) 2 and |Vr a( I 2 ( see Lemma 4.3). Indeed, 

. in5 a (-!) = {|0-f|<(l + ar)-*} 

14,1 W + |Vr Q(a±) ^| 2 = o( 1 + (ar)4+ 2* T ) 
is dominated by > 1 1+ °^ ) P 4+g ; 
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. in = fl >(1 + or)"*} 

l^| 2 (^) 2 + |Vr Q( . ±) ^| 2 = 0( r 



+ (ar)' 



is dominated by h' 2 a > \ a 1 +(° 8 r) ( 3 28) > fq^yi- 
Checking that v meets the conditions for a subsolution is absolutely analogous. In view of Lemma 
4.4, d Xg w > and similarly d xg v > 0. □ 

5. The solution. Existence and regularity. 

We have at our disposal a globally defined classical strict subsolution V to (1.1) lying below a 
classical strict supersolution W both of which are monotonically increasing in xg (in fact, strictly 
increasing in their interphases Kf n ). In this section we will explain why this engenders the existence 
of a classical solution u to (1.1), trapped in-between. Moreover, the solution will inherit some of the 
nice properties of the barriers V, W, such as monotonicity in xg and graph free boundaries F + (u) 
and F~(u). 

We will construct u as a global minimizer of I, constrained to lie between V and W. 

Definition 5.1. A function u £ Hl oc (JBL n ) is a global minimizer of I , constrained between V < W if 
for any bounded right cylinder 17 C R 9 

I(u, Q) < I(v, fi) for all v £ H 1 ^) such that V < v <W and u- v G H^(Q,). 

As usual, we obtain a global (constrained) minimizer u as a sequence of local (constrained) min- 
imizers on expanding bounded domains. For the purpose, we will verify that local minimizers are 
Lipschitz continuous with a universal bound on the local Lipschitz constant. This is done in the 
spirit of [10]. 

Afterwards, we will show that a global minimizer u which, in addition, meets certain simple 
geometric constraints, is actually a classical solution to our free boundary problem. This is achieved 
almost for free - by applying the regularity theory of minimizers to the energy functional Iq, developed 
in [10] and [11], to the functions 1 ± u. 

5.1. Existence of a local minimizer. Let ft C K.™ be a cylinder 

Crm = {x = (a/, x n ) e x K : \x'\ < R, \x n \ < h} 

and consider the minimization problem for the functional 

i(v,n)= f \vv\ 2 + i H <i, 

Jn 

where v ranges over the following closed convex subset of i? 1 (J7): 

S{Q) = {v e H\n) : V < v < W a.e.}. 
Let us show that there exists u € S(fl) for which the infimum of /(•, 17) over 5(17) is attained. 
Proposition 5.1 (Existence of monotone local minimizers). There exists u € S'(17) such that 

I(u,Q)=m:= inf 7(w,17). 

ves(n) 

Moreover, given that V and W are monotonically increasing in the x n -variable, u can also be taken 
to be monotonically increasing in the x n -variable. 
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Proof. For convenience use the simplified notation I(v) — I(v,Q). Obviously, the infimum m is 
non-negative and finite: 

< m < C := mxa(I(V),I(W)). 
Take a sequence Uk E S(fl) such that Cq > I(uk) \ m. Then 

IKIlHi = Klli 3 + l|Vu fc ||i 3 <|fi| + c , o 

is uniformly bounded, so by compactness we can extract a subsequence (call it Uk again) such that 

u/~ — > u in L? and a.e. and Vitfc — * Vu weakly in I? 

for some u £ S(Q). Claim that I(u) = m. It suffices to show that / is lower semicontinuous with 
respect to the weak-iJ 1 topology, i.e. 

I(u) < Uminf I{u k ) (5.1) 

which is done analogously as in [2]. 

We can produce a minimizer, which is monotonically increasing in the £ n -variable by applying 
a rearrangement. A monotone-increasing rearrangement in the a; n -variable, f —> f* satisfies the 
following properties (cf. [15]): 

(1) If / is monotonically increasing in the ir ra -variable, /* = /. 

(2) The functions / and /* are equimeasurable, i.e. |/ _1 (0)| = |(/*) _1 (0)| for any open interval 
OCR. 

(3) The mapping / — > /* is order-preserving, i.e. if / < g then /* < g* . 

(4) If / G H^CR.h), then /* G H^Crm) and 

iivrni 2 <iiv/n 2 L2 . 

Since V, W are monotonically increasing in the x n -variable, V* — V and W* = W; thus V < 
u* < W by order preservation under rearrangements. Moreover, u* G 1 (O) , so that u* G S(fl) and 
because of properties 2 and 4 above, 

m < I(u*) < I{u) = m. 

Thus, u* is a minimizer to / over S(tt), monotonically increasing in the £ n -variable. □ 

5.2. Lipschitz continuity of local minimizers. Employing standard arguments, we first establish 
continuity of local minimizers before we prove Lipschitz continuity with a universal bound on the 
local Lipschitz constant. 

We adapt the technique of harmonic replacements used by [2]. 

Definition 5.2. The harmonic replacement of u in the ball B C £1 is the unique function v G H 1 ^) 
that is harmonic in B and agrees with u on £l\B. 

Let 

b v ,w = Vin(v) u n in (w) 

and note that V is subharmonic in By.w whereas W is supcrharmonic in By.w- 

Below we show that the function u, constructed in Proposition 5.1, is continuous in By,w ■ 

Proposition 5.2 (Continuity). Let D d By.w ^= ^- Then the minimizer u, constructed in Proposi- 
tion 5.1, is in a Holder class C a (D) for some a > 0, depending on D. In particular, u is continuous 
in By,w ■ 
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Proof. Let Bp C D denote a ball of radius p, centered at some fixed point in D. Let v r be the harmonic 
replacement of u in the concentric B r C i? p . Since B r C Bywi where are subharmonic and 

superharmonic, respectively, the weak maximum principle implies that V < v r < W a.e. in B r . Thus 
v r € 5(f2) and /(u) < I(v r ). Therefore, 

/ \\7(u-v r )\ 2 = [ (|Vu| 2 - |Vu r | 2 ) < 2|B r | = c r n V0<r<p. 

J B r J By 

for some dimensional constant cq. Whence a standard dyadic argument in the spirit of [17, Theorem 
5.3.6] yields 

■f |Vm| 2 < C{1 + p- 1 ){\ + log 2 {p/r)) V0<r<p, 

J Br/4 

from which the statement of the proposition follows as in [17, Theorem 3.5.2]. 

□ 

Corollary 5.1. The function u is harmonic in f2j n (u) = {\u\ < 1}, subharmonic in {u < 1} and 
superharmonic in {u > —1}. 

Proof. From the previous proposition we know that u is continuous in By,w, therefore f2; n (u) C By.w 
is an open set. Thus for any x £ f2i n (w) we can find a small enough closed ball B = B r (x) C Sl m . Let 
v be the harmonic replacement of u in B. Since, \u\ < 1 on dB, the maximum principle implies that 
\v\ < 1. Combining the latter with the fact that harmonic extensions minimize the Dirichlet energy, 
we get that I(u,B) > I(v,B). However, by minimality, I(u,B) < I(v,B). Hence, I(u,B) = I(v,B), 
which in turn implies that 

l|Vn||| 2(B) = \\Vv\\l HB) 

So, u is itself the minimizer of the Dirichlet energy, meaning that u is harmonic in B. Since x £ Oj n 
is arbitrary, we conclude that u is harmonic in f2j n . 

The fact that u is subharmonic in {u < 1} and superharmonic in {u > —1} now follows from the 
mean- value characterization of sub/super-harmonic functions. □ 

Before we proceed to establish Lipschitz continuity, let us state the following definition related to 
the geometry of the pair of barriers V, W: 

Definition 5.3. We call the subsolution-supersolution pair (V, W) nicely intertwined in the 

bounded domain Q if 

F+{w)nn cn m (y) and F-(V)ntt cn m (w). 

so that F + (W)nfl stays a positive distance away from {V = ±l}nf2 and F~ (V)Hfl stays a positive 
distance away from {W = ±1} n ft. We say that V, W : R™ — > K are nicely intertwined globally 
if for every R > 0, there exists an ho — ho(R) large enough, such that (V,W) is nicely intertwined 
in all cylinders Q = Cn_h for h > ho(R). 

Proposition 5.3. Let D <s D' <s Q be compactly contained cylinders and suppose {V, W) is nicely in- 
tertwined in fl . Then there exists a constant K, depending onn, d{dD 1 dD'), d(F + (W)(1D' , F~(V)(~} 
D') and the Lipschitz constant of V,W in D' , such that \Vu\ < K in D. That is, u is Lipschitz- 
continuous in D. 

Proof. Since u G H 1 ^) and |u| < 1, 

Vu = Vuli u i<i a.e. 
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F+(V) 
F+(u) 
F+(W) 



F-(V) 
F-{u) 
F~{W) 



Figure 1. The free boundaries of a nicely intertwined pair (V, W) in a cylinder VL. 

Thus, it suffices to bound the gradient at points x <G Oi n (u) n D. Let B r = B r (x ) be the largest 
ball contained in i7; n n D' , centered at x$. We may restrict our attention to the situation when 
r < d ( dD ^ D ) ■ f or otherwise, using the gradient estimate for harmonic functions, 

\^0)\<^£n< 1( ^yy 

Obviously, in the situation when r < d ( dD ^ dD ) ; B r must touch the free boundary F(u) and not the 
fixed boundary dD' . 

Assume B r (xo) touches F + (u) at a point X\. We consider two cases determined by how close x\ 
is to F+(V). 

• Assume x\ is relatively close to F + (V): 

d{x!,F + (V)) = \x 2 - Xi\ < r/2 for some x 2 G F + (V). 

Note that \x 2 — xq\ < 3r/2 < d(D, D'), thus the segment between x$ and x 2 is contained in 
D' . The Lipschitz continuity of V in D' yields 

1 - «(a;o) < 1 - V(x ) < \\VV\\ L oo m \x2 - x \ < \\VV\\ L oo (D ,^r/2. 

Because 1 — u > is harmonic in B r {x$), Harnack's inequality implies that 

1 - u < c(l - u)(a;o) < c'||W||i,oo(rK)r in B r/2 . 

Hence, by the gradient estimate for harmonic functions we get the desired 

|Vu|0z o ) = |V(1 - u)\(x ) <-i (1 - u) < C\\VV\\ LO o {DI) . 

r JB r/2 (x„) 

• Assume that d{x 1 ,F+{V)) > r/2. Certainly, B r/2 { Xl ) C D', as r < d{dD,dD')/2. We may 
also assume that r < L = d(F + (W) n D', F~(V) n D'), for otherwise the gradient estimate 
for harmonic functions will immediately give us 

|Vu|(a;o) < j. 

With these assumptions in mind we see that B r / 2 (xi) C n in (V) n D' , so that u > — 1 on 
B r/2 (xi). 




32 



NIKOLA KAMBUROV 



Let v be the harmonic replacement of u in B r / 2 (xi). By the strong maximum principle 
\v\ < 1, so by minimality, 



/ \V(u-v)\ 2 <[ (l H<1 -l M<1 ) = |B r/2 (x 1 )n{u=l}|. 

JB r/2 (x 1 ) JB r/2 



(5.2) 



Now the argument for Lipschitz continuity of [2] goes through. It is based on the following 
bound for the measure of B r / 2 (x±) l~l {u = 1}: 

\B r/2 (x x )^{u = l}\(i (l-v)) 2 <Cr 2 [ |V(u-«)| 2 . (5.3) 

K JdB r/2 (x 1 ) ' JB T/2 {x 1 ) 

Hence, (5.2) and (5.3) imply 

f (1 - v) < Cr. 

JdB r/2 (x 1 ) 

Let 23 be a point on the segment between Xq and X\, which is at a distance r/4 from X\. 
According to Corollary 5.1, u is superharmonic in B r j 2 (x\), so 1 — u(x 3 ) < 1 — v(x 3 ). On 
the other hand, using a Poisson kernel estimate 



l-v(x 3 )<Cf (l-v)<Cr 

JdB r/2 (xi) 

for some dimensional constant C. Then by Harnack inequality, 

sup (1-u) < C'(l-u(x 3 )) < C'r. 

B4. r / 5 (x ) 

Applying a gradient estimate, we can conclude |Vu|(a;o) < G. 
The case when the ball B r (xo) touches F~(u) is treated analogously. □ 

5.3. Construction of a global minimizer. Take an increasing sequence of cylinders = Cn k: h k 
with Rk,hk 00 and let Uk be the minimizers to I(-,flk) over S(Qk) constructed in Proposition 
5.1. If (V, W) is a nicely intertwined pair, Proposition 5.3 implies that (for all large enough k) Uk are 
uniformly Lipschitz-continuous on compact subsets of 1". Therefore, one can extract a subsequence 
(call it again {uk}) which converges uniformly to a globally defined, locally Lipschitz continuous 
function u : M. n — > K, so that in addition 

Vu k Vu weakly in L~ c (M n ). 

Proposition 5.4. Assume that V,W is a pair of a globally defined subsolution and super solution 
to (1.1), which are monotonically increasing in x n and nicely intertwined with V < W. Then the 
locally Lipschitz-continuous function u : W l — > R, produced above, is monotonically increasing in the 
x n -variable, satisfies V < u < W and is harmonic in {u > 0}. Moreover, for any cylinder £1 C K™ u 
minimizes /(•, f2) among all competitors v G S(Cl) such that v — u £ H^(fl). 

Proof. The first three properties follow from the uniform convergence — ¥ u on compact sets. Let 
us concentrate on the minimization property: assume that there exists a cylinder O and a competitor 
v e S(Q), v - u e ffo(fi) such that 

I(v,ft) < I(u,n) - 5 
for some S > 0. Denote by Aft (CI) the t-thickening of fl: 

= :€ R™ : dist(s, Q) < t}. 
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For k large enough so that SI <s £lfe construct the following competitor v k ■ Cl k — > K for u k : 

fv(x) in Cl 

^ _ d^My i{x) + M Ufc (a;) in A t :=M(^)\n 
u fc (x) in ci k \Af t (ci) 

for some small enough £ > which will be chosen later. It is easy to check that v k £ S(Cl k ), therefore 

0<I(v k ,Cl k )-I(u k ,Cl k ) < 

< (l(v, Cl) - I{u, Cl)) + CI) - I(u k , CI)) + (l(v k ,A t ) - /(tt fc , A t )) 

<-6+ (I(u, CI) - I(u k ,Cl)) + (l(v k ,A t ) - I(u k ,A t )). (5.4) 

By the lower semicontinuity of /(•, CI) there exists a subsequence u kl such that 

I(u,Cl) - I(u kl ,Cl) <S/2. 

Now we claim that we can choose t so small that for all I large enough 

I(v h ,A t )-I(u kl ,A t )<5/2 (5.5) 

which will lead to a contradiction in (5.4). Indeed, |Vufe| < K is uniformly bounded on some fixed 
large cylinder CV D Aft (Ci) , so 

I(u kv A t ) < (K + l)\A t \ <Ct. 

Also, \Vu\ < K on Ct' 

d(x,Cl) 



|V Vfc; | 



V (u 



t 



'-(u kl - ufj 



< \Vu\ + ^\u kl -u\ + \Vu kl - Vu\ 



< 3K+j\u kl -u\. 



Thus, if e fe = sup n , \u k - u\ 

I(v h ,A t ) < C'\A t \{l + e k Jt+{e k Jt) 2 ) < Ct. 

for all I large enough, so that e k[ < t. Thus, if we choose t < S/(2C), the estimate (5.5) will be 
satisfied for all I large enough. □ 

5.4. Regularity of global minimizers. As mentioned in the introduction, there is an intimate 
connection between the energy functional / and the standard one-phase energy functional 

I («,fi) = / |V W | 2 + 1 {U>0} ueH'iCl). 
Jn 

as well as between the notions of viscosity (sub-/super-) solutions to (1.1) and (1.6). Recall, 

Definition 5.4. A viscosity solution to (1.6) is a non-negative continuous function u in CI such that 

• Au — in Cl p (u); 

• If there is a tangent ball B to F p (u) at some Xq £ F p (u) from either the positive or zero side, 
then 

u(x) — (x — Xq, v) + + o(\x — Xq\) as X — >• Xq, 

where v is the unit normal to dB at xq directed into Cl p (u) . 

Equivalently, a viscosity solution cannot be touched from above by a strict classical supersolution or 
from below by strict classical subsolution at a free boundary point. 

Definition 5.5. A viscosity subsolution (resp. supersolution) to (1.6) is a non-negative continuous 
function v in Cl such that 

• Am > in Cl p (u); 
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• If there is a tangent ball from the positive side B C £l p (v) (resp. zero side Q n (v)) to Fq(v) 
at some Xq £ F Q (v) and v denotes the unit inner (resp. outer) normal to dB at Xq, then 

v{x) = a(x — Xq, v) + + o(\x — x \) as x — > Xq, 
for some a > 1 (resp. a < 1). 
If a is strictly greater (resp. smaller) than 1, then v is called a strict viscosity subsolution (resp. 
supersolution). 

Remark 5.1. Indeed, suppose u minimizes I(£l) among all v £ H 1 ^), such that V < v < W , where 
V is a subsolution and W is a supersolution to (1.1) and let D £ be a (regular enough) domain 
such that either 

L>n-jV = -l} = or Dn{W = 1} = 0. 
In the first case, we readily see that 1—u minimizes Io(vq, D) among all admissible 1 — W < vq < 1 — V, 
with 1 — W being a subsolution and 1 — V - a supersolution to (1.6) in D. Similarly, in the second 
case, uq + 1 minimizes Iq(vq,D) among all admissible V+1<vq<W+1, with V + 1 being a 
subsolution and W + 1 - a supersolution to (1.6) in D. 

Below we will collect the regularity results concerning constrained minimizers of Iq, developed by 
[10]. For completeness, we will lay out the natural sequence of establishing regularity: starting from 
weaker notions and ending with the optimum, classical regularity. 

When a viscosity solution u to (1.6) arises from a minimization problem, u exhibits a non- 
degenerate behaviour at the free boundary in the sense that u grows linearly away from its free 
boundary in the positive phase. To make that statement precise we need the following definition. 

Definition 5.6. A continuous nonnegative function u is non- degenerate along its free boundary F p (u) 
in fi if for every G <e VL, there exists a constant K = K(G) > such that for every xq £ F p (u) fl G 
and every ball B r (xa) C G, 

sup u > Kr. 

B r (x ) 

On the way to establishing strong regularity properties for free boundary F p (u) one needs certain 
weaker, measure-theoretic notions of regularity. 

Definition 5.7. The free boundary F p (u) satisfies the density property (D) if for every G <s f2 there 
exists a constant c — c{G) > such that for every ball B r C G, centered at a free boundary point, 

\B r nn p (u)\ 

c - m - c - 

Here we should also recall the notion of nontangentially accessible (NTA) domains [14], which 
admit the application of the powerful boundary Harnack principles. 

Definition 5.8. A bounded domain D C K n is NTA if for some constants M > and ro > it 

satisfies the following three conditions: 

• (Corkscrew condition). For any x G dD, r < ro, there exists y — y r (x) € D such that 
M~~ 1 r < \y - x\ < r and dist(y,dD) > M~ l r; 

• (Density condition). The lebesgue density of D c = M. n \D at every point x € D c is uniformly 
bounded from below by some positive c > 

\B r (x)\ 

• (Harnack chain condition) If x\,x-i € D, dist{xi,dD) > e, i = 1,2 and \x\ — x 2 \ < me 
there exists a sequence of N = N{m) balls {B r .^ =l in D, such that x\ £ B ri , x 2 G B TN , 

< dist(B rj , dD) < Mr , j = 1, . . . , N. 
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We can now state the regularity results proved in [10] and [11] concerning constrained minimizers 
of I . We will say that the triple of functions (V , u , W ) defined on a vertical right cylinder SI — Cb„k 
satisfies the hypotheses H(51) if 

• Vq is a strict classical subsolution and Wq is a strict classical supersolution to (1.6) in f2 such 
that Vq < W , d Xn V > in {V > 0} and d Xn W > in {W > 0}. Moreover, F (V ), F (W ) 
are a positive distance away from the top and bottom sections of SI. 

• The function uo is monotonically increasing in x n and minimizes Iq{-, ^) among all competi- 
tors v € if 1 (SI) such that v - u £ H^Sl) and Vq < v < W in SI. 

Theorem 5.1 ([10], [11]). If (Vq,uq,Wq) satisfies the hypotheses H(fl), then: 

• [10], u is Lipschitz- continuous and non-degenerate along its free boundary F p (u); 

• [10], F p (u) satisfies the density property (D); 

• [10], F p (u) touches neither F p {Vq) nor F p (Wq); 

• [10], u is a viscosity solution to (1-6); 

• [10], For any vertical cylinder D <s SI the positive phase D l~l Sl p (u) is an NTA domain; 

• [10], The free boundary F p (u) n Cm h is given by the graph of a continuous function <p, 
F p (u) = {{x',x n ) : \x'\ < 3R/4,x n =(j>(x')}. 

• [11], // m&x\ x i\ <3 ji/4 \4>(x)\ < h — e, and e <C R < h then 

sup | V0| < C, 

|x'|<e/2 

where C depends on the dimension n, the Lipschitz constant of u, on h, e and the NTA 
constants of fl p (u) n C-m/ih- By the work of Caffarelli [(>], this implies 4>{x') is smooth in 
{\x'\ < e/4}. 

Let us revert our attention to the original problem. According to Proposition 4.4, we are in 
possession of a pair of a strict classical supersolution W : M 9 — > M and a strict classical subsolution 
V : R 9 — > E, such that V < W, both are monotonically increasing in the xg-variable (strictly 
increasing in that direction when away from their ±1 phases), and are, in addition, nicely intertwined 
(see Definition 5.3). By Proposition 5.4, we can then construct a globally defined, monotonically 
increasing in xg function u : K 9 — > K, such that u minimizes /(•, SI) among v G S'(Sl) for any vertical 
cylinder SI. Taking into account the observations we made in Remark 5.1, we can utilize the regularity 
results (Theorem 5.1) once we simply show that around every free boundary point x+ £ F + (u) there 
exists a vertical cylinder Sl + 3 x + , such that Sl + n {V = —1} = and around every free boundary 
point X- £ F~(u) there exists a cylinder Sl_ 3 X-, such that Sl_ n {W = 1} = 0. Furthermore, we'll 
need to show F + (u) n S7 + and F~ (u) n S!_ stay a positive distance away from the top and bottom 
of Sl_|_, respectively Sl_. This is the content of the next lemma. 

Lemma 5.1. Let V, W be the strict subsolution/ supersolution provided by Corollary 4-1 and u - the 
function constructed in §5.5. For every y' £ M 8 there exist an y n £ R and an R > small enough 
such that 

-\ < V(x',y n ) < W{x\y n ) < ~ 
for all x' £ R 8 with \x' -y'\<R. Thus, if 

V+iv') = {(x',x n ) : \x' -y'\ < R, < x n - y n < h} 
tt-(y') = {(x',x n ) : \x' -y'\ < R, < yn-Xn < h} 

and h > is large enough, the monotonicity of V and W in the x n -direction guarantees that 

• CL+(y') n {V = -1} = and F+{V),F+{W) exit from the side o/Sl+; 

• Sl_(y') n {W = 1} = and F-{V),F~{W) exit from the side o/S7_. 
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In particular, (1 — W, 1 — u, 1 — V)(x', —x n ) satisfy hypotheses H(il + (y / )) , while (V + 1, u + 1, W + 
l)(x',x n ) satisfy hypotheses H(Q,_(y')). 

Proof. Fix y' G E 8 . Since W(y',x n ) = ±1 for all large positive (negative) x n , there certainly exists a 
y n such that W(y',y n ) = \. For a small enough R > 0, we can ensure < W(x',y n ) < § whenever 
|x' — y' | < i?. Since, VF and V were constructed so that 

< W{x) - V(x) < - VieR 9 (Corollary 4.1) 

we see that 

-\<V{x' lVn )<W{x\x n ) < l - \x'-y'\<R. 

□ 

Localizing at any free boundary point, we immediately invoke the lemma above and the regularity 
results (Theorem 5.1), establishing the desired Theorem 1.1. 

Appendix A. Supersolutions for Jr^ 

The two objectives of this appendix are 

• To state the results of Del Pino, Kowalczyk and Wei [12] concerning supersolutions for the 
linearized mean curvature operator on i*^: 

V(j> (VFoo-V^VFocA 



H'[F 00 ]{cj>) = ± 



H[F 00 +t ( f>}=di v 

t=o 



yi + iv^p (i+ivfooi a ) 3 / 2 ; 

where H[] is the mean curvature operator (MCO), and to describe the relation between 
H'[Foo] and the Jacobi operator Jp^, on Too. 
• To obtain the refined estimate 3 = O ( ^ . 

We recall that in polar coordinates (r, 9), Foo(r,6) = r 3 g(0) and the function g{9) is smooth and 
satisfies: 

(1) g(6) = -,g(f - 6) 



near 



(2) s (0) = (0-f)(l + ps(6l -!)* + ...) 

(3) g'(6)>0 for(9e(0, §) and g'(0) = s'(tt/2) = 0. 
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A.l. The relation between the Jacobi operator and the linearized MCO. For a domain 
U C M 8 , let S : U -> K be an arbitrary C 2 (U) function and denote by E = {(#', x 9 ) e [/ x K : x 9 = 
its graph. Denote the standard projection onto K 9 by 

7T : M 8 x R ->■ M 8 . 

We can identify functions 4> defined on U with functions on E in the usual way: 

0S = 0° 7T. 

We'll abuse notation and use the same symbol cf> to denote both. A long, but straightforward 
computation yields the following interesting formula relating the linearized MCO associated with S 
to the Jacobi operator Js := As + \ Az\ 2 on the graph E: 

■fc ( , ) - mm) - v ™>- vs , . (a.i) 

Note that if E is a minimal graph, i.e. i? [5] = 0, we recover the well-known relation 

J s ( . I = if'[S](<£). (A.2) 
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For our purposes, we would like to estimate the size of the error term in (A.l) when S = F^. 
A.2. Computation of H[Foo] and iV^-Foc])!- First, we compute iffi^]. Since div (jff 2 ^) = 0, 



ff^oo] = div 



-div 



VF a 
\VF n 



-div ^ 



IVFool^l + IVFooPdV^ool + V 71 + 1^00 | 2 ) 



• V 



i \ = v-Foo vg 



where Q(x') := y/l + \VF^(\VF ao \ + ^/T+W F ^\ 2 ~)■ Note that Q is bounded from below by 

Q{x') > 2\VF QO \ 2 = 2r 4 (9.g(0) 2 +g'{9) 2 ) > 2mr\ 
where m — xmng e [ 0t „/2] (9g(#) 2 + d'{9) 2 ) > 0- Also, 



|V|V^|| 2 = (2r7vTF) 2 + (r9'^==) 



Iv^i + IVFooi 2 ! 2 = 



|VF a 



VlVF, 



v/1 + lVFo.l 2 



= |VF„ 



0(r 2 ) 

: |V|Vfoo|| 2 
1 + IV.Fool 2 " 



0(r 2 ). 



Thus, 
and 



|VQ| -0(r 3 ) 



IVQI 

Q 2 



| J ff[ J F 1 o ]|<^=0(r- b ). 



To compute |V-ff[-Foo]|i observe that 



The last two summands are obviously 0(r~ 6 ). Let us bound 
IVCVFoo • VQ)| 2 - (^((foo)rQr + r~ 2 (F^eQe)) % 

+ r- 2 (d e ((FocWr + r- 2 ^)^^)' 



Because of (A. 3), Q r ,r 1 Qg = 0(r 3 ). Furthermore 



Foorri T F^q^q^t F&Q , 



0(r), 



and 



Q rr ,r Q r9l r Qee = 0(r ). 



Thus, |V(VFoo • VQ)| 2 = (3(r 8 ) and the first summand in (A. 4) is then 



V(VFoo-VQ) „/r 4 



= O 



.10 



as well. We conclude 



iVi^lQ 2 

iV^FooD^OCr" 6 ) 



= 0(r- 6 ), 



(A.3) 



(A.4) 



(A.5) 
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A. 3. Supersolutions for Jr x ■ In [12, §7.2] the authors study the linearized MCO H'lFao] and show 
that it admits two types of supersolutions away from the origin. We will call those Type 1 and Type 
2 in parallel with the labels we used in Section §4.2. Because of (A. 5), formula (A.l) becomes 

so that one can then cook up supersolutions for the Jacobi operator Jr^ . 

• Type 1 supersolution for the linearized MCO (cf. the proof of Lemma 7.2 in [12]): 

There exists a smooth function <f>i = cj>\(r,9) = r~ e qi(9) with qi(6) > and even about 
9 = 7r/4 that satisfies the differential inequality 

tf'^oo] (</>!) < r>r 
for sufficiently large ro- Thus, (A.6) implies that hi — — ?== = € C^fT^) satisfies 

Jr^h 1 <--^ + 0(r-*-<) < - rT i iT7 (A.7) 

for sufficiently large r > tq. 

• Type 2 supersolution for the linearized MCO (cf. the proof of Lemma 7.3 in [12]): 

For every | < r < | there exists a function </> 2 = fair, 9) — rq 2 (9), defined for 8 € {7r/4 < 
9 < tt/2}, such that 

tf'^cc]^) < ^S(|,|], r>r 

for sufficiently large ro. Moreover, q 2 (9) is smooth in (f , ^] and has the following expansion 
near 9 = j + : 

qi{6) = (9- ^) T (a + a 2 (9- j) 2 + ■ ■ ■ ) where a Q > 0. 
Therefore, /it = , ^ 2 = satisfies the differential inequality 



1 

in € (f , f ] for sufficiently large r > r$. 



^ 2 <-^+0(r-W)<-^ (A. 



A. 4. Computation of -ffoo,3- We will compute the second fundamental form A of the graph Too 
and then estimate the sizes of the principal curvatures. 

Let y = (ur cos 9, vr sin 9, r 3 g (9)) € r m , with i,ieS 3 cR 4 and consider local parametrizations 
u(tt, ti, ts), and v(si, s 2 , S3) of S 3 around u and v, respectively, such that 

u(0)=u d u u(0) = Ti, i= 1,2,3 
w(0)=w d Si v(0)=ai, i = 1,2,3 

where are orthonormal bases for T^S" 3 and T V S 3 , respectively. Then 

P(r,9,ti,Si) = (ur cos 9,ir sin 9,r 3 g(9)) (A.9) 

defines a local parametrization of near y. 

In the system of coordinates {r, 0,i,,s,} the metric tensor near y takes the form g = g2 © g S ym, 
where 

l + 9rV 3r 5 53 ' \ _/ (r 2 cos 2 9)U 3 



g2 ~' 3r 5 . 93 ' r 2 (l+ 5 'V) J ,gsym ~ ' (r 2 sin 2 9)V 3 >' {XM) 
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In the expression for g sym above, U^jVs are 3x3 matrices that depend only on {tk},{sk} with 
^3(0) = ^3(0) = ^3, the identity 3x3 matrix. We will also need the inverse of g, g _1 = g^ 1 © g^ n , 
where 

„4„/2 o„3„ n > 



-1/ s 1 / 1 + rV 2 -3rW A ifn 2 , 2 . 

12 f (cos 9)-*I 3 \ 
m(V) { (sin0)- 2 / 3 )' 



(A.11) 



The unit normal is given by 

v(y) = — i=((Foo) u u, (Foo) v v, -l) = 

V 17 

= — -= (r 2 (35 cos 0-5' sin r 2 (3.g cos 6 + g' sin 0)5, -l) . 

We calculate the second fundamental form (A^ij — —diP ■ djV — dijP ■ v at y to be Aoo — 
(^00)2 © (^oo)sym, where 

/ . \ 1 \ 1 ( 6rg 2r 2 g' 
^^IW r3(3 3 + 3 ") 

M s f s _ J_ / rcos0(F (X3 )„7 3 

\At V rsm ^(Foo )„J 3 

The principal curvatures of Too at ?/ are the eigenvalues of the matrix Aoog -1 , i.e. the eigenvalues 
Ml, ^ 2 (each of multiplicity 3) of (A 00 ) sym (g sym )~ 1 : 

M2 = ^^ = ^(3 5 + «/cot0) 
V<7 r sin ft y (j 

and the eigenvalues Ai, A 2 of 

1 ( 0{r{9-D) 



(A 00 ) 2 {g 2 )- 1 = 



a 3 / 2 V * 4 0(r 5 (6»- f)) 



4 - 



Since g' cotO = 0(1) and g'tanO = 0(1), we see that /^i and ^2 are 0((1 + r) ). Note further that 

Mi + M2 = -£= (63 + 2g' cot 20) = 0( V ) ■ ( A - 12 ) 

On the other hand, Ai, 2 = 0((1 + r) _1 ) as well, since 

Ai + A 2 = Trace((A 00 ) 2 (g 2 )- 1 ) = o(^^) (A.13) 

AiA 2 =det(A 00 ) 2 det(g 2 )- 1 -0(1)^- =0(r- 6 ). 

Now (A.12) and (A.13), combined with the fact that the principal curvatures are all of order 
Oiil + r)- 1 ) imply that 

Hoo,3 = 3(mi + M2)(M? - M1M2 + Ma) + ( A i + A 2)( A ? - A i A 2 + Af.) = o( ^r^ ) • 
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Appendix B. Bounds on the gradient and hessian of h(r,6) on r^. 

Once again we will make use of the {r, 9, ij, Si} system of coordinates (A. 9) in order to estimate 
the first and second covariant derivatives of a function h = h(r, 9) € C 2 (r oo ) which depends only on 
r and 9. 

Lemma B.l. The gradient and hessian of h = h(r,9) € C 2 (r oo ) satisfy: 

\D roo h\ = 0(\d r h\ + {r~ l d + r- a )\d e h\) (B.l) 
\D^h\ = o(\8 2 h\ + (r" 1 !? + r-3)\d 2 r0 h\ + (r~H + r-' i ) 2 \djh\^ + 

+ r- x 0(\d r h\ + {r~ l ti + r- 3 )\d g h\^ + o((r~ 6 + r~ 2 tf 2 ) 
where $ := \8 — 7r/4| and the constants in the O-notation depend on g. 



deh 



(B.2) 



Proof. In order to carry out the computations, we adopt the standard Einstein index notation. That 
way, we write 

\D T ^h\ 2 = h l h, and \D 2 Ji\ 2 = hjh 3 \ 
where i,j range over the list of coordinates {r, 9, {s/c}} 

hi = dih, h l =g lk h k with g« = (g-%- 

and 

hi =dihP +T\ k h k . 

ik 



In the expression above, V 3 -, are, of course, the Christophel symbols: 



r « = (di&k + dkgu - d lglk ) . (B.3) 

Since h = h(r,9), we have {Dr^hl 2 = h r h r + h hg. Using (A. 11) we calculate 

1 + rV 2 3r 3 gg' 
n = n r tig 

h = -^h r+ r ^ +9 ^ 2 h e 
a a 

and taking into account that g(9) = 0(d) and 9g 2 + g' 2 is uniformly bounded from above and from 
below by positive constants, we conclude 

\Dr^h\ 2 = 0(\h r \ 2 + r-H\h r \\h e \ + (r~ 3 + r^d) 2 ^ 2 ) 

so that (B.l) is verified. 

The computation of the hessian is slightly more involved. For convenience we will denote by Greek 
letters a, ft, 7, etc. indices that correspond to coordinates r, 9, and by Latin Z,m,n, etc. indices that 
correspond to coordinates {U},{si}. First, note that the "cross term" contribution hjhf 1 = 0, 
because 

hi a = d t h a + r« bP = 0, 

as dih a = and 

rg, = \g ai {d lglP + d pgh - d lgw ) = 0. 

Thus, iD^h] 2 = S 2 + S syia , where 

S2 := ha^hp", S syra := hi m h rn l . 



Let us first deal with S'sym. We see that 



°sym — 1 l a L m(3 11 a ! 
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where 

1 la = \Olg ma + O a gl m - O m gl a ) — -g d a 

Noting that 

9 a (r 2 sin^)/3 ; ' 

we obtain 

A more refined estimation of h e gives 



, 5 a (r 2 cos 2 6»)/ 3 

9 Q gsym(2/) = fl /J .J, 



S sym = (^) 2 ^ + 2/i r //-(cot(9-tan(?) + (/i e ) 2 3(cot 2 (9 + tan 2 0). 



ft 9 = o 

and since 



/ |(|ftr|/r + (r- e + r- a l? a )|/ lfl /s / |)) 
cot 9 - tan0 = 0(1/ g') cot 2 + tan 2 6 = 0(1/ g' 2 ) 



we conclude 



S sym = 0(\h r \ 2 r- 2 + r~H 2 \h e \ 2 + (r- 6 + r^ 2 ) 2 \h e /g'\ 2 ). (B.4) 



sym 

Now we proceed with the computation of 5*2. For the purpose we need to calculate the Christophel 

7 

a/3- 



symbols T'la- The derivatives of g2 are 



_/ 36rV 15r 4 ff5 ' \ _/ 18r 4 ff5 ' 3r 5 (g' 2 + gg") 

drg2 ~ { 15r 4 55 ' 2r + 6 ff 'V ) 6 ° g2 \ 3r 5 ( 5 ' 2 + gg") 2r*gg> 

which we then plug in (B.3) to obtain: 

K r \ = g^_ / d rgrr \ = 6 / 3r 3 . 9 2 \ = / 0(r~H 2 ) 

T e rr J 2 \ 2d r g er - deg rr J ° \ r 2 gg' J \ 0( r - 2 $) 

T r 9r \ = &2^[ degrr \ 1 ( 6rV \ ( 0(0) 

Tl ) 2\d r g ee ) a \ + 3r 3 (g' 2 + 3 ff 2 ) J \ 0(r^) 

T ee \ _ &2 1 ( ld e g re - d r gee \ = 1 / r 5 (3«? ff " - g' 2 ) -r\ = ( 0(r) 

degee ) v\ r 4 g'(g" - 3g) ) \ 0(d) 



For convenience define the following expressions that measure the magnitude of the first and second 
derivatives of h: 

J 7 {h a ):=\h r \ + ('dr- 1 + r- 3 )\he\ 
S(d 2 ap h, h a ) := \d 2 r h\ + (fir- 1 + r~ 3 )\d 2 r h\ + (tfr" 1 + r- 3 ) 2 \d 2 h\ + ^^(K) 

A straightforward computation yields: 

h r r = d r h r + r rr h r + r r r8 h e = o{s(dl P h, ha)) 

h e e = d e h 9 + T e er h r + T r ee h e = 0(S(d 2 a/3 h, h a )) 

h r e = d r h e + r e rr h r + r%h = 

= 0(r- J hg + r- & d 2 r9 h + tir^Sidlph, ha)) + 0($r- 2 F(h a )) = 
^0((r- 3 + ^r- 1 )S(d 2 a ^h, h a )) 

h e r = d e h r + r gr h r + r gg h 6 = 

= O(d 2 r0 h + r-^d 2 h + §h r + r~%) + 0(dh r + (r^ 5 + ^V" 1 )^) = 
= 0((r- 3 + dr- 1 )- 1 S(d 2 ap h 1 ha)), 
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whence we conclude 




(B.5) 



□ 
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